AD-A044  801 
UNCLASSIFIED 

SPERRY  RESEARCH  CENTER  SUDBURY  MASS 
IMPULSE  RESPONSE  TARGET  STUDY* (U) 

SEP  77  C L BENNETT*  R HIERONYMUS*  H MIERAS  F30602-76-C 

RADC-TR-77-273 

F/6  17/9 

-0209 

NL 

1 1 cf2 

1 A0A04480( 

B 

• 

u 

^ - -- 

i 

1 

1 

y 

i 

t 

1 

1 

<1 

— 

j 

i 

) 

I 

1 

1 

i 

1 

\ 

1 

\ 

' 

J 

1 

K 

1 

\ 

^ 

1 

1 

■ - 

4 

p 


Q 

IMPULSE  RESPONSE  TARGET  STUDY 


Sperry  Research  Center 


Approved  for  public  release;  distribution  unlimited. 


I >" 

Q- 

I O 

! O 


• ~ — 

O u. 


ROAAE  AIR  DEVELOPMENT  CENTER 

Air  Force  Systems  Command 

Griffiss  Air  Force  Base,  New  York  13441 


D D O ^ 


I 


This  report  has  been  reviewed  by  the  RADC  Information  Office  (01)  and 
is  releasable  to  the  National  Technical  Information  Service  (NTIS) . At  NTIS 
it  will  be  releasable  to  the  general  public  including  foreign  nations. 

This  report  has  been  reviewed  and  is  approved  for  publication. 


APPROVED 


C/jOSEPH  L.  RYERSON 
Technical  Director 
Surveillance  Division 


FOR  THE  COMMANDER: 

JOHN  P.  HUSS 

Acting  Chief,  Plans  Office 


If  your  address  has  changed  or  if  you  wish  to  be  removed  from  the  RADC 
mailing  list,  or  if  the  addressee  is  no  longer  employed  by  your  organization, 
please  notify  RADC  (DAP)  Griff iss  AFB  NY  13441.  This  will  assist  us  in  main- 
taining a current  mailing  list. 


Do  not  return  this  copy.  Retain  or  destroy. 


REPORT  DOCUMENTATION  PAGE 


2 GOVT  ACCESSION  NO 


1 REPORT  number 


R/\DC-TR-77-27  3 


4 TlTUF  (tuid  Subtitle) 


IMPULSE  RESPONSE  TARGET  STUDY 


7 AUTHORf*; 

C.  L.  /flennett, 
K.  Hieronymus 
11 . . Mleras 


» PERFORMING  ORGANIZATION  NAMF  AND  ADDRESS 


Sperry  Research  Center 
100  North  Road 
Sudbury  >U  0177b 


1 CONTROLUNG  OFFICE  name  AnD  ADDRESS 

Rome  Air  Development  Center  (OCTN) 
Criffiss  APR  NY  134A1 


14  monitoring  agency  N AMF  A ADOR^S^;/  dtlfrtont  froui  Cux»Uoil»rnJ  OUtce) 


Same 


16-  DISTRIBUTION  STATEMENT  (ol  fhi.i  Keporl) 


Approved  for  public  release;  distribution  unlimited 


17  DISTRIBUTION  STATEMENT  lof  thf‘  aSsRA.  I cnR-raU  iri  filocli  JO,  It  ihtlvtfnt  froni  Hfpott) 


Same 


8 supplementary  notes 
RADC  Project  Ensineer 


Daniel  L.  Tauroney  (0C1>I) 


reverse  side  il  necessarv  and  ld<*nfllv  bv  hlorfc  number) 


19  KEY  WORDS  (Coiittni 


Space-Time  Integral  Equation;  Impulse  Response  Augmentation  Technique 
Target  Ceometry;  Impulse  Response;  Frequency  Response 


10  abstract  (ConUniie  «>/i  reverse  side  it  itecenserv  and  identitv  bv  Mo<  A nfnbar' 

The  space-time  integral  equation  approach  is  extended  to  the  problen  ol  sin-pie 
aircraft  shapes.  The  target  consisted  of  a cvlindrical  fuselage,  and  rectan- 
gular wings,  stahilixers,  and  rudder.  Calculations  were  compared  to  measure- 
ments. The  inverse  scattering  problem  was  ref or-uil ated  to  l ik.  advant.ige  of 
direct  time  domain  method  of  solution.  Solutions  for  the  splu>re,  sphere  capped 
cylinder,  and  flat  end  cylinder  were  developed. 


EDITION  OF  I NOV  fill  18  OHSOLE  'F 


RKAD  INSTRUCTIONS 
BK.hOKK  COMPBKTINC.  KOKM 

3 BFClPltNT'S  catalog  NUMBER 

5 Type  of  K£Hort  a pchiod  coy 

RED 

Final  Technical  Report. 
15  Apr  76  - 39  Apr  77. 

OR9:  RFPORT  number  I 

N/A 

B contract  or  grant  NUMBERfsl 

F30602-76-C-0209-  ' ' 

10  program  El  ement  project  , task  I 

AREA  & WORK  UNIT  NUMBERS 

62702F 

45060575 

12  REPORT  D_^ae 

September  L977 

NUMBER  OF  pages 

122 

15  security  class,  fol  fhi*  reporfp 

UNCLASSIFIED 

ISa  DF-CL  ASSIFIC  ATION  DOWN  GRADING 

[n/a 

- 

TAHU-:  ')!■  onti-:nt:' 

Soc^tion  Page 

] INTRODUCTION  1 

2 SPACE-TIME  INTEGRAL  EQUATION  SOLUTION  EOR  SIMPLE  AIRCRAFT 

MODELS  3 

2.1  General  Problem  3 

2.2  Space-Time  Integral  Equation  Solution  5 

2.2.1  Derivation  of  Space-Time  Integral  Equation  5 

2.3  Numerical  solution  12 

2. ‘1  Time  Domain  Moa.surements  17 

2.4.1  Time  Domain  Meaauremont  System  19 

2.4.2  Comparison  of  Measured  and  Calculated 

Responses  21 

2.5  Calculated  Response  of  Aircraft  Model  21 

3 EXTENSION  OF  IMPULSE  RESPONSE  AUGMENTATION  TECHNIQUE  31 

3.1  Review  of  the  Impulse  Response  Augmentation  Technique  31 

3.2  Application  of  Cylinder  with  Fins  37 

3.  3 Results  ‘54 

4 TIME  DOMAIN  INVERSE  SCATTERING  - AXIAL  INCIDENCE  55 

4.1  Derivation  of  Equations  56 

4.2  Numerical  Solution  59 

4.3  Results  64 

5 TIME  DOMAIN  INVERSE  SCATTERING  - OBLIQUE  INCIDENCE  73 

5.1  General  Problem  73 

5.2  Derivatroti  of  Inverse  .scattering  Equations  74 

5.3  First  Order  Ramp  Resi'onse  of  Curvilinear  Patch  B1 

5.-1  Numerical  Solution  Procedure  B7 

5.5  Inversi'  -Scatteiinq  Results  B8 

6 SUMMARV  AND  CONCLUSIONS  95 

7 REFERENCES  '^9 

8 APPENDIX  50.1 

B.l  Improved  Space-Time  Integral  Equation  Solution  Techni>.jue  101 

8.2  Useful  Di  fferiuit  i.il  tieomet  ry  Relat  ionships  108 

8.3  Integration  fiver  S'(’l  f-Pat  cli  HO 

8.4  Transformat icm  of  Self-l'atc!i  Correction  to  the  Incident 

H Coordinate  .System  114 

8.5  Projected  .Area  of  Nivvily  Illumin.ited  Regicni  llL 

8.6  Initial  Ramp  Rer.ponsu'  of  .Specular  Point  1 -'0 

I i i , ■ 1 v 


LIST  OF  ILLUSTRATIONS 


Figure  Page 

1 Aircraft  sketch.  4 

2 Simple  aircraft  model.  4 

3 Equivalent  problem  of  scatterhng  from  surfaces  with  fins  attached.  6 

4 Geometry  of  far  scattered  field.  18 

5 Targets  used  for  comparison  of  calculated  and  measured  responses,  22 

6 Comparison  of  calculated  and  measured  responses,  an  = 1,  a = 90°.  23 

7 Comparison  of  calculated  and  measured  responses,  an  = 0.5, 

a =90°.  24 

8 Comparison  of  calculated  and  measured  responses,  a = 0,  an  = 0.5.  25 

9 Calculated  response  of  aircraft  model  at  0°,  30°,  60°,  90°.  26 

10  Calculated  response  of  aircraft  at  90°,  120°,  150°,  and  180°.  27 

1 1 Comparison  of  responses  of  several  finned  cylinders  at  axial 

incidence.  29 

12  Impulse  response  augmentation  technique.  35 

13  Geometry  of  cylinder  with  fin.  38 

14  Smoothed  impulse  response  of  finned  cylinder  at  various  angles 

of  incidence  for  TE  polarization.  ’ 39 

15  Smoothed  impulse  response  of  finned  cylinder  at  various  angles 

of  incidence  for  TE  polarization.  40 

16  TE  response  of  finned  cylinder  with  radius  a for  0°  angle  of 

incidence.  45 


17 


TE  response  of  finned  cylinder  with  radius  a for  30°  angle  of 


incidence. 


r ' 


46 


TO 


-J  ^ C/5 


(5 


I 

~~ — ->v  j 


V 


LIST  OF  ILLUSTRATIONS  (Com.) 


Figure  Page 

18  TE  response  of  finned  cylinder  with  radius  a for  60'^  angle 

of  incidence  48 

19  TE  res()onse  of  finned  cylinder  witfyradius  a for  90“  angle 

of  incidence  ' 49 

20  TE  response  of  finned  cylinder  with  radius  a for  120“  angle 

of  incidence.  51 

21  TE  response  of  finned  cylinder  with  radius  a for  150°  angle 

of  incidence.  52 

22  TE  response  of  finned  cylinder  with  radius  a for  180°  angle 

of  incidence.  53 

23  Geometry  of  rotationally  symmetric  scattering  problem.  60 

24  D irect  time  domain  solution  procedure.  62 

25  Geometry  of  objects  used  for  test  of  inversion  procedure.  65 

26  Far  field  ramp  response  (backscatter  direction)  for  four  objects 

tested.  66 

27  Coordinate  axes  and  view  angle  for  perspective  plots.  67 

28  Contour  estimeates  for  a sphere.  68 

29  Contour  estimates  for  a cylinder  with  two  sphere  caps.  69 

30  Contour  estimates  for  a cylinder  with  a flat  front  end  and  a 

sphere-capped  back  end.  70 

31  Contour  estimates  for  a cylinder  with  a sphere-capired  front  end 

and  a flat  back  end.  71 

32  Inverse  scattering  of  rotationally  symmetric  surface.  75 

33  Unknown  region  A S for  body  with  rotationaly  symmetry, 

to()view.  63 


VI 


LIST  OF  ILLUSTRATIONS  (Com.) 

Figure  Page 

34  Self- patch  integration  at  point  inside  A S,  83 

35  Inverse  scattering  results  — cross-section  of  sphere  quadrant.  89 

36  Inverse  scattering  results  — sphere.  89 

37  Inverse  scattering  results  — cross-section  of  prolate  spheroid.  91 

38  Inverse  scattering  results  — prolate  spheroid.  92 

39  Inverse  scattering  results  — cross-section  of  sphere-capped 

cylinder.  93 

40  Inverse  scattering  results  — sphere  capped  cylinder.  94 

41  Projected  area  of  newly  illuminated  region.  117 


VII 


EVAI.UATION 


The  contractor  was  successful  in  computing  the  impulse  resijonsc-  for  a 
simple  aircraft  target  at  several  aspect  anciles.  The  si)ace-time  integral 
approach  has  led  to  manv  new  rerults  associated  with  this  model.  The.-t-  re- 
sults can  be  ust;d  to  obtain  the  radar  response  of  those  targets  over  the 
entire  spectrum.  This,  in  turn,  can  be  used  to  compute  the  response  of  a 
taruot  due  to  any  incident  radar  waveform,  regardless  of  wave  shape  or  carrier 
freauency.  The  technique  developed  on  the  effort  not  only  applies  to  th.e 
aircraft  identification  portion  of  TPO  R2E  but  also  to  the  analysis  of  tar 
gets  which  are  constructed  of  composite  materials  and  have  radar  absorbing 
materials  incorporated  in  highly  reflecting  areas. 
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SECTION  1 


INTRODUCTION 

This  document  is  siibmitted  as  the  final  report  in  response  to  tiie 
1‘1,'qui  rements  set  forth  in  Contract  No.  F30602-76-C-0209  between  the  Sperry 
Coiporate  Research  Center  (SCRC) , Sudbury,  Massachusetts,  and  tlie  Air  Force 
Systems  Command,  Rome  Air  Development  Center  (RADC) , Griffiss  Air  Force  Base, 
New  York.  The  primary  objectives  of  this  proqram,  performed  during  the 
period  15  April  197G  through  30  April  1977  were  to  (1)  improve  and  expand 
computer  programs  tliat  had  previously  been  delivered  to  RADC  and  (2)  to  ex- 
tend the  techniques  previously  developed  by  SCRC  for  computing  tlie  smoothed 
impulse  response  and  the  impulse  response  of  relatively  simple  targets  to 
complex  targets  which  have  tins,  edges,  and  curved  surfaces.  The  first  ob- 
jective was  carried  out  and  described  in  an  interim  report  [1].  The  second 
objective  is  discussed  in  detail  in  this  report. 

In  Sec.  2 the  space-time  integral  equation  approach  is  extended 
to  the  case  of  simple  aircraft  models.  Techniques  are  developed  for  co;npu- 
ting  the  smoothed  impulse  response  of  a simple  aircraft  model  consisting  of 
a cylindrical  fuselage,  rectangular  wings,  rectangular  stabilizers,  and  a 
rectangular  rudder.  The  techniques  are  verified  by  comparison  with  direct 
time  domain  measurements.  Smoothed  impulse  response  results  are  displayed 
for  the  simple  aircraft  model  at  numerous  aspect  angles  and  specific  portions 
of  the  response  arc  related  to  scattering  by  local  regions  of  the  target. 

This  appears  to  be  the  first  solution  technique  that  has  been  developed  to 
liandle  tl.ose  sim^ile  aircraft  models. 

I’tie  extension  of  the  impulse  response  augmentation  technique  to 

cylinders  with  fins  attached  is  discussed  in  Sec.  3.  In  jjarticular,  results 

are  obtained  for  the  example  of  a cylinder  with  square  fins  and  a length  to 

o o 

diameter  ratio  of  2:1  tor  a series  of  aspect  angles  that  go  from  0 to  180  . 

In  Sec.  4 the  inverse  scattering  problem  for  a.xial  incidence  on 
a rotational ly  symmetric  target  is  reformulated  to  take  advantage  of  tlie 
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iitoct  time  domain  method  of  solution.  The  previously  developed  iterative 
method  was  cumbersome  and  time  consuminy.  The  new  "marchintj  on  in  time"  ap- 
[■roach  is  faster  and  more  straightforward.  In  Sec.  4 this  new  approach  is 
described  and  demonstrated  on  a sphere,  sphere-capped  cylinder  and  a sphere- 
cap,  flat-end  cylinder  and  displayed  on  three  dimensional  plots.  The  results 
are  compared  with  those  obtained  using  the  iterative  technique. 

The  extension  of  the  direct  time  domain  approach  to  the  inverse 
scattering  problem  to  the  case  of  rotational ly  symmetric  targets  at  oblique 
incidence  is  discussed  in  Sec.  5.  The  techniques  are  developed  and  demon- 
strated on  three  targets.  Conclusions  are  presented  in  Sec.  6. 
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SECTION  2 


[ 


SPACE-TIME  INTEGRAL  EQUATION  SOLUTION  FOR  SIMPLE  AIRCRAFT  MODELS 

The  space-time  integral  equation  is  valid  for  any  conducting  target 
shape  and  solution  techniques  have  been  developed  for  a number  of  important 
target  classes  (2-Gl . In  this  section  the  solution  techniques  are  developed 
for  a simple  aircraft  model  and  verified  by  comparison  of  calculated  and 
measured  results.  Sec.  2.1  describes  the  general  problem  and  Sec.  2.2  des- 
cribes the  space-time  integral  equation  and  its  solution  for  this  new  target 
class.  The  time  domain  measurement  procedure  is  reviewed  and  the  measured 
and  calculated  results  are  compared  in  Sec.  2.3.  The  calculated  smoothed 
impulse  response  of  the  simple  aircraft  model  is  presented  and  discussed  in 
Sec.  2.4. 

2.1  GENERAL  PROBLEM 

The  target  class  considered  in  this  work,  consists  of  a conducting 
surface  with  protruding  wings  as  shown  in  Fig.  1.  The  problem  is  to  cal- 
culate the  smoothed  impulse  response  of  this  target  class  as  a function  of 
aspect  angle.  In  order  to  make  the  problem  amenable  to  solution  by  existing 
space-time  integral  equation  techniques,  the  target  is  modeled  as  shown  in 
Fig.  2.  In  this  model  the  fuselage  is  represented  by  a conducting  right- 
circular  cylinder  and  the  wings,  stabilizers,  and  rudder  are  represented 
by  thin,  flat  rectangular  conducting  plates. 


In  this  work  a smoothed  impulse  excitation  is  treated,  yielding 
the  target  response  well  into  the  resonance  region.  High  frequency  augmen- 
tation techniques  exist  to  extend  the  solution  over  tlie  entire  frequency 
domaiti.  These  techniques  have  been  described  earlier  (4J  and  in  Sec.  3 of 
this  report. 

In  the  numerical  implementation  described  here,  tl\e  direction  of 
incidence  is  taken  to  bo  in  the  plane  of  symmetrv'  aiivl  tlie  ]iolar ization  is 
taken  to  be  perpendicular  to  the  plane  of  sv'mmetry  (TF.  mode)  as  illustrated 
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in  Fig.  2.  It  should  be  noted  that  most  of  these  r<  t !■ 'ti-  ir»-  t in- 
herent to  the  space-time  integral  equation  technique  and  ^ i.t-  r ived  at 

the  expense  of  only  computation  time  and  complexity.  For  exaniif,  t:.<.  fuse- 
lage could  be  any  convex  surface,  the  various  fins  could  be  flat  i lates  with 
more  general  straight-edge  contours,  and  the  symmetry  condition.  ss)ull  be 
removed  to  allow  computation  at  arbitrary  incidence  with  arbitrary  jolariza- 
tion. 

2.2  SPACE-TIME  INTEGRAL  EQUATION  SOLUTION 

The  problem  of  determining  the  scattering  by  cylinders  with  fins 
attached  is  of  great  practical  interest,  since  this  serves  as  a model  for 
numerous  missiles  and  aircraft.  To  the  authors'  knowledge,  there  has  been  no 
adequate  procedure  developed  for  the  solution  of  this  complex  scattering 
problem  in  either  the  time  domain  or  the  frequency  domain.  This  secti  -.i 
describes  the  method  for  solving  the  problem  of  scattering  by  a cylinder 
with  multiple  sets  of  fins.  This  is  an  extension  of  techniques  developed 
earlier  for  cylinders  with  a single  set  of  fins  [5] . Basically  this  approach 
consists  of  developing  two  simultaneous  space-time  integrodif ferential  equa- 
tions and  their  subsequent  computer  solution  by  marching  on  in  time.  These 
equations  contain  terms  which  may  be  interpreted  as: 

(a)  The  influence  of  cylinder  currents  on  other  cylinder 
currents . 

(b)  The  influence  of  fin  currents  on  cylinder  currents. 

(c)  The  influence  of  fin  currents  on  other  fin  currents. 

(d)  The  influence  of  cylinder  currents  on  fin  currents. 

The  neighborhood  along  the  line  where  a fin  is  attached  to  the  cylinder  is 
accounted  for  by  application  of  boundary  conditions  at  the  edge  of  the  fin. 

2.2.1  Derivation  of  Space-Time  Integral  Equation 

The  technique  used  here  to  develop  the  space-time  integrodif feren- 
tial equation  for  surfaces  witli  fins  attached  is  to  consider  tlie  equivalent 
of  this  scattering  problem  shown  in  Fig.  3.  In  this  equivalent  statement 
of  the  2^toblem,  the  conducting  surface  has  been  replaced  by  source  n.urface 


currents  and  the  conducting  fins  have  been  replaced  by  source  fin  currents. 

Since  these  source  currents  are  now  radiating  in  free  space,  the  techniques 

used  for  the  case  of  scattering  from  surfaces  alone  can  be  combined  with  the 

techniques  used  for  the  case  of  scattering  from  flat  plates  alone  to  solve 

the  problem  of  scattering  from  surfaces  with  fins  attached.  The  integral 

equation  for  the  surface  currents  is  obtained  by  considering  the  problem  in 

•> 

Fig.  3,  in  which  the  scatterer  has  been  replaced  by  source  currents  J which 
are  equal  to 

. -»■ 

J = a X 11  . 
n 


Since  these  currents  are  radiating  in  free  space,  the  free-space  Green's 
function  may  be  used  to  find  the  expression  for  the  total  H field  at  an  ar- 
bitrary point  in  space,  which  is 


H(r , t) 


-vi  ->- 

H (r,t) 


j'  _1 3 

4tt  .1  ^2  Rc 


J(r ',!) 


X a ds ' 
R 


T = t - R/c 


(1) 


whore 

r = position  vector  to  the 
observation  point 

► ^ 

r = position  vector  to  the 
integration  point 


-►  -»■ 

r - r 


c = the  speed  of  light. 


An  integral  equation  for  the  current  density  J may  be  obtained  by  speci.iliz- 

ing  the  arbitrary  space  p'oint  r to  a point  on  the  surface  of  the  scatterer 

-►  *• 

and  then  applying  the  boundary  conditions  to  cast  II  in  terms  of  J.  Fer- 
forming  the  limiting  procedure  and  applying  the  boundary  conditon  yields  the 
integral  cejuation 
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J(r,t)  = 2a  X H (r,t)  + 
n 


2tt  J % -^^^2  Rc  3i  ' ' 


T)  X a MS' 
R 


T = t - R/c  ( 2 ) 


where  r is  now  located  on  the  surface  of  the  scatterer.  Equation  (2)  is 
a space-time  vector  integral  equation  for  the  current  density  J on  the  sur- 
face of  the  scatterer.  The  first  term  in  the  right-hand  side  of  Eq.  (2) 

may  be  considered  the  source  term  and  represents  the  direct  influence  of  the 

>• 

incident  field  on  the  current  at  the  observation  point  (r,t).  The  integral 
term  on  the  right-hand  side  of  Eq.  (2)  gives  the  influence  of  currents  at 
other  surface  points  on  the  current  at  (r,t). 

The  factor  of  2 in  going  from  Eq.  (1)  in  free  space  to  Eq.  (2)  on 

the  surface  arises  due  to  the  following  well  known  boundary  value  argument. 

Consider  the  small  patch  at  the  surface  point  r.  The  field  H at  r is  due 

to  the  effect  of  the  patch  plus  the  effect  from  the  rest  of  the  surface  and 

-*• 

the  incident  field.  If  we  move  r just  inside  the  conducting  surface  (where 
we  know  the  field  to  bo  zero),  the  contribution  from  the  patch  changes  sign. 
The  magnitude  of  the  effect  of  the  patch  is  thus  the  same  as  that  due  to 
t)\e  remainder  of  the  surface  and  the  incident  field.  Taking  the  limit  as 
the  patch  size  goes  to  zero,  we  thus  obtain  for  the  field  just  outside  the 
surface  twice  the  value  of  Eq.  (1) . Performing  the  cross  product  with  a^ 
we  then  obtain  Eq.  (2) . 

The  currents  flowing  on  the  surface  are  a function  not  only  of  tlic 
incident  field  but  also  of  the  currents  at  other  cylinder  surface  points  and 
fin  surface  points.  Equation  (2)  can  be  expanded  to  display  these  two  con- 
tributions explicitly  and  give  the  total  space-time  integral  equation  for 
the  cylinder  surface  currents  as 


. 2a__ « f . -i  ♦ i sr  ■>=' 

R ! 

C 


+ J (r  ,T)  X dS 

2TT  ' n 1,2  R 3t  F R 


T = t - R 


T = t - R 


r 


In  Iaj.  (3)  the  units  of  time  have  been  chanqed  from  seconds  to  light-meters, 
where  one  light-meter  is  defined  as  the  time  it  takes  an  electromagnetic  wave 
traveling  at  velocity  c to  travel  a ilistance  of  one  meter. 


The  integrodi t ferential  equation  for  the  fin  currents  is  obtained 
by  starting  with  tlie  expression  for  tlie  electric  field  in  terms  of  the  elec- 
tric and  magnetic  potentials.  The  geometry  of  the  problem  is  displayed  in 
Fig.  2.  In  tliis  particular  formulation,  the  fins  lie  in  the  x-z  and  y-z 
planes  and  the  incident  field  makes  an  angle  a with  the  x-z  plane.  For  tliis 
geometry,  the  total  electric  field  is  given  by 


>■  >•  ->-i  ->• 

E(r,t)  = E (r,t) 


( r , t ) 

^ ~Tt — 


V4)(r,t) 


(4) 


where 


A = magnetic  vector  potential  such  that  H = V ^ A 
= electric  potential 
M = permeability  of  sjjace. 

Next,  apply  to  liq.  (4)  tlte  Lorentz  gauge  relation 


V • A + i 4^  = 0 

I’t 


where  r is  the  permittivity  of  space,  yielding 


+ V (V  • A)  - tit- 


3t 


3t 


Dt‘ 


(5) 


Tlieti  applying  tlu'  E-field  boumiary  conditioTi  tliat  tlie  component  of  the  total 
E-field  tangent  to  the  fin  vanisher.  everywhere  on  the  fin.  Eg.  (A)  can  be 
expii's.st'd  as 


\,  V • A ( r , t ) 


i 

l_Mr,t)  'E  (r,t) 
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or,  in  a more  simplified  form 


2-*-  -*■ 

n A +fiiA 


= -e 


at 


(7) 


where 


A ^ + 


rv  ^ A 

Jx  dt  / 

the  wave  operator. 
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T = t - R 


In  order  to  obtain  an  expression  for  the  fin  currents  directly,  it 
is  necessary  to  express  the  magnetic  vector  potential  A as  the  sum  of  two 
components,  one  due  to  the  observation  point  ("self"  patch)  and  the  other 
due  to  the  remainder  of  the  patches.  If  the  fin  grid  patches  are  small 

enough,  then  the  current  over  a aiven  oatch  can  be  assumed  constant.  "’his 
makes  it  possible  to  express  A as 


A(r,t)  = Y Jp(r,t)  + — 


; j(r  ,T) 
V R 


Non-Self 

Patches 


dS 

t - R 


(8) 


where 


y = 2i 


As 


, TT 


6 ' 

2J 


6 

2/' 


As  = area  of  patch  containing  r 
6 = thickness  of  fin. 


Substituting  Eq.  (8)  into  Eq.  (7)  yields  the  fin  currents  directly 
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where 
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f 

J R 

Non-Seif 

Patches  t = 


ds  ” 
t - R 


The  surface  currents  contributing  to  are  both  cyiinder  currents  and  fin 

currents.  The  contribution  of  both  of  these  can  bo  separated  in  Eg.  (9) 
and  dispiayed  expiicitiy  to  give  the  expression  for  the  current  flowing  on 
tlie  f in 


2 ► >•  .•  + 1 I 9e^ 

’ (r,t)  +-X!  J (r,t)  = - - c 

f l-  1 dt 


9 > > V 

A^,  (r,t)  + A (r,t) 

r C 

NS 


- TC  A^  ( r , t ) + A ( r , t ) 

F i- 

NS 


(10) 


wliere  A^(r,t)  is  the  magnetic  vector  potential  due  to  tlie  cylinder  current. 


The  fin  currents  at  the  free-space  edges  and  at  the  fin-cylinder 
join  are  given  by  the  boundarv'  conditions 


wliere 


J (r,t)  = 0 , J„(r,t)  " (free-space  edge) 

(r,t) 

— i— = 0 , J =0  (fin-cylinder  join)  (11) 

Art 


J,  (r,t) 


current  component  {perpendicular  to  tlie 
edge  or  join 


(r,t) 


current  component  {parallel  to  the  edge  or 
join. 
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Thus,  tlie  space-time  integral  equations  (3)  and  (10)  together  with  the  boun- 
dary conditions  in  Eq.  (11)  form  the  solution  for  the  surface  currents  on 
a cylinder  with  fins  attached. 


In  most  applications,  the  quantity  of  interest  is  the  far-scattered 

field.  This  can  be  computed  directly  from  the  surface  currents.  Tlie  far- 

Tt^s 


field  expression  for  the  scattered  magnetic  field,  H , may  be  obtained  from 
Eq.  (1)  by  noting  that  t)ie  contribution  of  the  first  term  within  the  integral 
becomes  negligible  when  R becomes  very  large.  In  addition. 


and 


R r 


Thus,  the  expression  for  the  far  field,  with  time  expressed  in  light-meters, 
becomes 


->s  -> 

rH  (r,t) 


9d(r',T) 


47:  .1 


3l 


X a dS " 
r 


(12) 


T = t - R 


2.3  NUMERICAL  SOLUTION 


To  obtain  a solution  of  the  space-time  integral  equations  for  a 
cylinder  with  fins  attached,  each  of  tlie  components  of  the  problem  must  be 
represented  mathematically.  First,  the  geometry  of  the  scatterer  and  tlie 
characteristics  of  the  incident  field  must  be  numerically  specified.  From 
these,  the  surface  and  fin  currents  can  bo  computed  using  numerical  repre- 
sentations of  Eqs.  (3)  and  (10) . Tliis  is  accomplished  by  carrying  out  the 
integration  and  dif ferentiation  numerically  and  using  a "marching  on  in  time" 
procedure.  The  far-scattered  field  can  be  computed  directly  from  the  current 
densities  by  using  a numerical  representation  of  Eq.  (12). 
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To  describe  tlie  scattercr  geometrically,  the  scattering  surface 
is  divided  into  curvilinear  patches  of  approximately  equal  area  with  a space 
sample  point  at  the  center  of  eacli  patch.  The  spacing  of  these  sample  points 
(and  thus,  the  size  of  the  patches)  is  chosen  small  enough  to  give  both  a 
good  representation  of  the  scatteror  itself  and  of  the  currents  that  exist 
on  the  scatterer.  The  sample  point  spacing  affects  the  time  increment  At 
at  whicti  the  current  densities  can  be  calculated.  The  time  increment  must 
not  exceed  tlie  time  it  takes  a wave,  moving  at  tlie  speed  of  light,  to  travel 
between  tlie  closest  space  points.  This  insures  that  the  integral  equation 
can  be  expressed  as  a recurrence  relation  in  time  and  that  a matrix  inversion 
is  not  necessary  to  obtain  a solution. 


The  si’ace-time  integral  equations  (3)  and  (10)  represent,  in  prin- 
ciple, the  solution  of  the  scattering  problem  for  cylinders  with  fins  at- 
tached for  an  arbitrary  incident  field.  Although  the  equations  can  be  solved 
for  each  incident  field  separately,  it  is  very  inefficient  to  do  so.  Since 
in  most  practical  scattering  problems  the  excitation  is  a plane  wave,  a more 
efficient  way  to  approach  this  problem  is  to  compute  the  scattered  response 
when  the  incident  wave  is  an  electromagnetic  impulse.  Once  the  impulse  re- 
sponse of  a target  has  been  obtained,  the  response  due  to  any  incident  plane 
wave,  whose  spectrum  is  contained  within  the  spectrum  of  an  impulse,  can  be 
calculated  by  a simple  convolution  procedure.  Moreover,  t)ie  impulse  resixjnsc 
is  intimatelv'  related  to  the  actual  geometry  of  tlie  target,  and  thus,  the 
potential  for  developing  tecliniques  to  determine  the  impulse  response  of 
a scatterer  by  an  inspection  of  its  geometry  is  ever  present. 


For  the  numerical  solution,  however,  it  is  not  practical  to  use 
an  ideal  impulse  for  an  excitation.  Thus,  in  this  solution  a regularized 
or  smoothed  impulse  is  used.  The  form  of  this  illumination  at  the  origin  is 
the  Gaussian  regularization  of  an  impulse,  namely 

llNt)  = (13) 

/it 

which  converges  to  a delta  functional  as  n goes  to  infinity.  The  time  do- 
main integral  equations  can  be  solved  exactly  for  bodies  with  linear 


A 


-13 


dimensions  up  to  several  pulse  widths  of  this  regularized  impulse.  In  this 
solutioti  consideration  is  limited  to  bodies  of  this  size. 

The  currents  flowing  on  the  cylinder  and  on  the  fin  are  computed 
using  numerical  representations  of  Eqs.  (3)  and  (10) . For  the  purpose  of 
discussing  the  numerical  solution,  only  the  x-components  of  these  equations 
will  be  considered.  Similar  representations  are  used  for  the  calculation  of 
tlie  other  components. 

For  the  numerical  solution,  Eq.  (3)  for  the  x-component  of  the 
cylinder  current  at  sample  point  i and  time  t is  represented  as 


^ 2 n - n ^ E F(J  ) n . + n . 


J (r . ,t)  = 

Cx  1 yi  z zi  y 2tt 

£=i 


n n 

_F 

'x'  ’’yi  R zi  R 


F(J  ) n . - F(J  ) n . As 

y yi  R z zi  R x. 


(14) 


where 


Y = self-term  correction  factor  for  observation  patch  i 

i 


(h\h^ 

y z 


N 

F(J) 


J 

X 


(r^,T) 


the  y-  and  z-components  of  the  incident  field  at  patch  i, 
time  t 

number  of  grid  patches  on  the  cylinder  and  fin 
J(rj^,T)/R  + 3j(r^,T)/ax 

the  x-component  of  the  current  density  at  patch  £,  time  T 


T = t - R 

t = time  in  liglit-meters 

R = distance  from  the  integration  patch  £ to  the  observation 

2 2 2 
patch  i:  R = /(x.  - X.)  + (y.  - y.)  + (z.  - z.) 

1 C 1 > 1 X 

a . = unit  normal  at  patch  i,  (n  . , n . ,n  , ) 
nr  XI  yi  zi 
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unit  vector  from  patch  i to  patch  i,  (n  , ii  ,n^  ) 

Rx  Ky  Rz 

X - x„ 

i y 

R 


R 

^ ~ 

R 

area  of  patch  i . 

The  time  differentiation  and  interpolation  necessary  for  the  eval- 
uation of  the  integrands  appearing  in  Eg.  (3)  are  performed  numerically  by 
representing  the  sux'face  current  with  a fourth-order  polynomial.  In  order 
to  achieve  the  best  accuracy,  the  five  points  used  for  tlio  representation 
arc  chosen  such  that  the  current  is  evaluated  as  near  as  possible  to  the  mid- 
dle of  them. 

Tlie  numerical  representation  of  Eg.  (10)  for  the  x-component  of 
tlic  fin  current  at  (r.,t)  is  obtained  by  representing  the  second  time  deri- 
vative of  the  fin  current  by  a three-point  difference  approximation,  yielding 
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At  = solution  time 
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X 111 

X 


step 


+ A 


NS 


-15- 


N 

)i=i 


J 

R 


ASj 


X = t - At  - R . 


X'hu  space  and  time  derivatives  in  Eq.  (15)  are  evaluated  using 
three-  and  fiv'e-point  difference  approximations.  For  the  time  interpolation 
a linear  app>roximation  is  used. 


The  equations  for  the  surface  and  fin  currents  are  solved  with  a 
digital  computer  by  simply  marching  on  in  time.  The  computation  starts  at 
a point  in  time  before  the  incident  field  reaches  the  scatterer  and  proceeds 
sequentially  in  time  in  the  same  manner  that  nature  would  solve  the  problem 
in  the  real  world.  It  is  imp>ortant  to  note  that  since  the  minimum  spacing 
between  space  sample  points  on  the  surface  is  not  less  than  At,  then  Eqs. 

(14)  and  (15)  give  the  current  density  in  terms  of  other  currents  at  times 
not  later  than  (t  - At)  wtiich  are  already  icnown.  Thus,  the  integral  Eqs. 

(3)  and  (10)  have  been  reduced  to  recurrence  formulas  in  time  and  the  need 
to  perform  matrix  inversions  has  been  eliminated. 


Once  the  current  densities  have  been  computed,  the  far-scattered 
field  can  be  calculated  directly  using  a numerical  representation  of  Eq. 
(12).  For  computation  purposes,  Eq.  (12)  is  expanded  in  rectangular  coor- 
dinates yielding 
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where 


a + a + a 

XX  y y z z 

a=n  a+n  a+n  a. 
r rx  X ry  y rz  z 

The  time  differentiation  and  interpolation  necessary  for  the  evaluation  of 
3j(r^,T)/9T  are  performed  numerically  using  a five-p>oint  Lagrange  interpola- 
tion formula. 


A sketch  of  the  geometry  parameters  used  for  the  description  of 
the  far-scattered  field  is  drawn  in  Fig.  4.  In  this  sketch  the  incident 
field  is  shown  proi^agating  in  t>ie  yz-plane  and  making  an  angle  (x  with  the 
z-axis.  The  scattered  fields  produced  by  the  surface  currents  are  then  com- 
puted in  the  two  principal  planes.  The  scattered  field  is  computed  in  the 

yz-plane  at  angles  ^ with  respect  to  the  direction  of  propagation  of  the 

yz 

incident  wave.  In  this  plane  the  two  orthogonal  components  used  to  represent 


the  scattered  field  are  the  component  perpendicular  to  the  yz-plane,  H 


-^s 


yzx 


and  tlie  component  tangent  to  tlie  yz-plane,  The  other  plane  in  v;hich 

the  scattered  field  is  computed  is  the  px-plane,  which  is  formed  by  tlie 

direction  of  propagation  of  the  incident  wave  and  the  x-axis.  In  this  plane 

the  two  components  used  to  represent  the  scattered  field  are  tlie  component 

s 

perpendicular  to  the  px-plane,  H , and  the  component  tangent  to  the  px- 
s ^ 

plane,  H . These  scattered  fields  are  computed  at  angles  d’  , which  are 
pxt  px 

measured  with  respect  to  the  direction  of  propagation  of  the  incident  wave. 


2.4  TIKF.  DOMAIN  MEASUREMENTS 

The  computer  program  incorporating  tiie  methods  described  in  the 
preceding  section  was  verified  by  comparing  the  calculated  time  domain 
smoothed  impulse  response  with  measurements  made  on  the  Sperry  Research  Cen- 
ter time  domain  scattering  range.  ^ This  patented  measurement  system  [7], 
whicli  has  been  described  in  detail  [5,8,9],  is  reviewed  in  Sec.  2.4.1  and 
the  calculated  and  measuredjK  esult.s  arc  compared  in  Sec.  2.4.2. 
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2.4.1  Time  Domain  Measurement  System 

The  time  domain  scattering  range  is  a system  for  obtaining  a low- 
noise  response  in  the  nanosecond  region.  It  consists  of  a ground  plane  scat- 
tering range,  a step  function  source,  a sampling  oscilloscope  receiver,  and 
a labc-atory  instrumentation  computer  for  control  and  processing.  The  system 
signal  source  is  a high-voltage  switch  which  generates  a 300  V step  function 
with  a risetime  less  than  100  ps.  The  signal  is  radiated,  virtually  undis- 
torted, from  a wire  transmitting  antenna  protruding  through  a circular  ground 
plane  20  foot  in  diameter.  This  wave  illuminates  the  target  and  the  result- 
ing scattered  waveform  is  received  on  a coaxial  horn  antenna,  which  essen- 
tially smootlies  and  differentiates  the  signal  and  thus  provides  the  smoothed 
impulse  response  of  the  target.  The  received  waveform  is  sampled  by  a 12 
GHz  oscilloscope  that  has  been  triggered  by  the  initial  pulse  and  whose  sam- 
pling gate  deflection  is  under  the  control  of  a small  instrumentation  compu- 
ter. Unprocessed  data  are  displayed  on  the  oscilloscope  CRT  while  the  sam- 
pled-and-held  waveform  is  passed  through  a low-pass  filter,  digitized,  read 
into  the  computer,  and  stored  on  magnetic  tape  automatically.  This  system 
has  been  designed  to  correct  the  long-term  timing  drift  and/or  amplifier 
drift.  In  addition,  the  waveforms  are  stored  in  such  a way  that  they  are 
ready  for  the  subsequent  operations  of  averaging  (to  remove  short-term  noise) 
and  baseline  processing.  The  effects  of  a time  varying  baseline  are  sub- 
tracted from  measured  waveforms  to  improve  system  accuracy. 

Tlie  salient  characteristics  of  the  range  are  the  speed  and  sim- 
plicity with  which  multi -octave  frequency-domain  data  can  bo  obtained.  Tliose 
advantages  accrue  because  Uie  time-domain  scattering  range  yields  an  "uncon- 
taminated" interval  of  time  between  the  arrival  of  the  direct  wave  ami  ttie 
arrival  of  unwanted  reflections.  Targets  are  usually  located  anywliere  from 
two  to  five  foot  from  the  transmitting  antenna.  Tlie  response  from  tlu'  anten- 
na tip  and  the  taiale  edge  occur  at  approximately  15  ns.  Thus,  a "clean 
window"  exists  between  4 ns  and  15  ns  which  can  be  used  to  view  the  target 
responses.  The  entire  region  between  tiie  direct  transmission  and  the  tal>le 
edge  respotisc  forms  a convonit'nt  t ime  "window"  to  view  the  target  response 
an>'i  allows  one  to  "gate  out"  (in  tliru')  unwanted  reflections.  Thus,  undis- 
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torted  transient  target  responses  can  be  viewed  without  resorting  to  elaborate 
and  expensive  anechoic  chambers.  In  addition,  a single  time-domain  measure- 
ment obviates  the  requirement  for  tedious  measurement  of  the  amplitude  and 
phase  responses  at  many  frequencies. 

The  accuracy  of  the  measurement  system  has  been  estimated  for  the 
results  presented  in  this  section.  The  peak  of  the  incident  pulse  as  meas- 
ured on  the  sampling  oscilloscope  is  approximately  400  mV,  and  a typical  tar- 
get response  has  a peak  value  in  the  vicinity  of  10  mV.  When  using  the  10  mV 
scale  on  the  sampling  oscilloscope,  the  standard  deviation  of  the  Scimple  mean 
is  estimated  to  be 

0-  = 0.5  mV 
v 

if  16  scans  are  averaged.  Thus,  the  estimated  standard  deviation  of  the 
sample  mean  v is  in  the  vicinity  of  5%  of  the  peak  value  of  the  target  re- 
sponse . 

The  signal  received  directly  from  the  source  is  also  measured. 

Taking  the  distances  between  source,  receiver  and  target  into  account,  the 
frequency  response  of  the  target  is  obtained  by  numerically  deconvolving  the 
source  signal  from  the  response.  The  source  is  not  exactly  the  Gaussian 
smoothed  impulse  used  for  the  calculated  response  (Eq.  13).  Therefore,  a 
convolution  procedure  is  carried  out  on  the  measured  response  to  permit  com- 
parison with  the  calculated  response  in  the  time  domain.  The  width  of  the 
smoothed  impulse  is  chosen  to  bo  slightly  greater  than  the  width  of  the 
source  signal.  This  has  the  effect  of  reducing  high  frequency  noise  in  the 
response.  Taking  the  smoothed  impulse  width  greater  yet,  i.e.  reducing  the 
higher  frequency  contents  of  the  incident  pulse,  smoothes  the  response  fur- 
ther at  the  expense  of  a reduction  in  the  fast  structure  of  the  response. 

In  these  measurements  the  incident  pulse  width  was  about  0.59  ns  or  7".  The 
comparisons  were  made  at  smoothed  impulse  widths  of  8"  and  16".  The  target 
lengths  were  8"  and  12". 
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2. “1.2  Comparison  of  Measured  and  Calculated  Responses 


The  targets  chosen  for  comparison  were  8"  long  and  12"  long  right 
circular  cylinders,  each  of  radius  2",  and  each  with  a set  of  4"  square  fins 
at  the  ends,  as  illustrated  in  Fig.  5. 

In  Fig.  6 the  responses  of  these  two  targets  are  shown  at  broadside 
o 

incidence  (a  = 90  ) for  an  = 1 . The  larget  radius  a is  2"  and  the  nominal 

pulse  width  is  defined  as  4/n  = 8".  The  time  scale  is  t/a  (t  is  in  light- 

meters)  and  the  response  has  been  normalized  to  r fi/a.  (See  ref.  4 for  a 

o 

discussion  on  scaling  and  normalization.)  The  solid  curve  represents  the 
calculated  response,  the  dotted  curve  the  convolved  measured  resixinse. 

While  the  agreement  is  remarl<.ably  close,  it  is  interesting  to  ob- 
serve the  deviations  (of  less  than  10%  of  pea)<  value)  in  the  rapidly  fluc- 
tuating peaces  of  the  signals.  These  disagreements  are  greatly  reduced  when 
the  comparison  is  made  at  an  = 0.5  (corresponding  to  a pulse  width  of  16") 
as  is  done  in  Fig.  7 for  these  same  targets.  Inaccuracies  thus  begin  to 
appear  at  shorter  wavelengths.  It  is  believed  that  these  inaccuracies  are 
almost  entirely  due  to  the  relatively  coarse  patch  structure  of  tlie  numerical 
representation  of  the  target  used  in  the  response  calculations.  On  the  basis 
of  experience  with  space-time  integral  equation  calculations  on  other  targets 
with  edges,  it  is  )cnown  that  even  greater  accuracy  can  be  obtained  by  more 
detailed  structuring  of  the  target  representation,  in  particular,  near  the 
fin  edges. 

As  an  additional  comparison,  tlie  results  are  given  in  Fig.  8 for 
tlic  same  targets  viewed  tail-on  (Ot  = 0)  with  an  = 0.5. 

2.5  CALCULATED  RESPONSE  OF  AIRCRAFT  MODEL 

Tlie  model  chosen  for  those  calculations  is  shown  in  Fig.  2.  The 
calculated  responses  are  plotted  in  Figs.  9 and  10  for  aspect  angles  from  0^ 
to  18o”  at  3o”  intervals. 

It  is  instructive  to  note  the  composition  of  the  resj-'Oiise.  Consiii- 


-21- 


FIG.  8 Comparison  of  calculated  and  measured  responses,  tv  = 0,  an  = 0.5. 
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er  the  tail-ori  aspect  (ni  = 0°)  : Tlie  initial  £>ositive  swing  is  due  to  t)ie 

return  from  the  cylinder  face  and  the  leading  edges  of  the  small  fins.  The 
negative  swing  at  t = -4  is  liue  both  to  the  differentiating  effect  of  the 
cylinder  face  and  the  trailing  edge  of  the  small  fins.  A creeping  wave  ef- 
fect due  to  small  fin  currents  interactions  via  the  cylinder  surface  can  be 
seen  near  t = -1.5  just  before  the  large  positive  peak  from  the  leading  ed- 
ges of  the  large  fins  at  t = 0.  This  is  followed  by  a large  negative  swing 
from  the  largo  fin  far  edges  and  fin-fin  creeping  wave  interactions  at 
t = 8.  The  creeping  wave  effects  continue  after  that,  damping  out  near 
t = 16.  In  Fig.  11  these  effects  are  illustrated  in  more  detail.  Shown 
are  the  responses  of  the  some  target  at  0°  aspect  without  fins  and  with  each 
set  of  fins  separately.  For  comparison  purposes,  the  time  scale  (in  light- 
meters)  of  tile  grapli  and  the  sketch  of  the  incident  pulse  and  the  space 
scale  of  the  target  are  the  same;  similarly  tlie  amplitude  of  incident  and 
response  pulses  are  the  same  scale.  The  behavior  described  above  is  very 
evident  in  this  series  of  graphs. 

Returning  to  Figs.  9 and  10  it  is  noted  tliat  as  the  aspect  changes, 

the  response  from  the  cylinder  face  moves  toward  t = 0;  while  the  response 

from  the  small  fins  moves  more  rapidly  toward  t = 0,  joining  that  of  t)ie 

o o 

large  fins  at  Ct  = 90  . At  a = 90  the  positive  peak  (due  to  the  fins)  near 
t = 0 is  strong  enough  to  obscure  the  negative  swing  expected  at  tliat  time 
from  a broadside  response  of  the  cylinder  alone.  The  negative  swings  of 
cylinder  body,  small  and  large  fins  combine  in  tlie  large  response  near 
t = +1.  After  t = 2 the  response  is  entirely  due  to  creeping  waves. 

At  angles  beyond  90°,  the  responses  from  the  cylinder  face  and 
the  large  fins  become  evident  before  t = 0 as  they  separate  with  increasing 
U,  while  the  response  from  the  small  fins  shows  as  a small  perturbation  at 
times  increasingly  later  tlian  t = 0. 
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SECTION  3 


EXTENSION  OF  IMPULSE  RESPONSE  AUGMENTATION  TECHNIQUE 

The  space-time  integral  equation  approach  solves  the  scattering 
problem  directly  in  the  time  domain.  This  approach  is  valid  for  any  excita- 
tion; however,  the  most  useful  excitation  has  been  found  to  be  the  regular- 
ized (or  smoothed)  impulse  given  in  Eq.  (13) . The  response  due  to  this 

s 

excitation,  r H it/a),  is  the  regularized  (or  smootlied)  impulse  response  of 
o f 

the  target  and  is  computed  exactly  with  the  space-time  integral  equation 
using  the  techniques  described  in  Sec.  2.2.  This  approach  yields  results 
for  targets  with  sizes  up  to  several  pulse  widths  or,  equivalently  in  the 
frequency  domain,  several  wavelengths. 

A technique  has  been  developed  which  uses  the  space-time  integral 
equation  approach  as  the  basis  and  extends  the  results  to  obtain  the  impulse 
response  and  frequency  response  of  an  arbitrary  target  over  the  entire  spec- 
trum [4] . This  technique  has  been  demonstrated  for  several  smooth  convex 
targets  including  the  sphere,  the  prolate  spheroid,  and  the  sphere-capped 
cylinder.  It  has  also  been  applied  to  targets  with  edges  including  right 
circular  cylinders,  flat-end  sphere-cap  cylinders  at  axial  incidence  and 
square  flat  plates  [4,5] . The  purpose  of  the  sections  which  follow  is  to 
describe  the  extension  of  the  impulse  response  augmentation  technique.  Sec. 
3.1  provides  a review  of  the  impulse  response  augmentation  technique.  The 
extension  to  cylinders  with  fins  attached  at  arbitrary  incidence  is  described 
in  Sec.  3.2. 

3.1  REVIEW  OF  THE  IMPULSE  RESPONSE  AUGMENTATION  TECHNIQUE 

In  order  to  simplify  the  notation  in  these  sections  the  electro- 
magnetic field  variables  are  equated  to  their  linear  system  counterparts  as 
follows : 

H^(t/a)  ->■  e(t)  = incidence  pulse 
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where 


# 


r (t  /a)  >■  r (t) 

o f 


h(t) 


t/a  ->•  t 


smoothed  impulse  response 

impulse  response 
time 


H^(ka)/a  ->■  E(oj)  = transform  of  e(t) 

r H^(ka)/a  R(a>)  = transform  of  r(t) 
o 

H(aj)  = frequency  response 
ka  0)  = frequency 


= incident  magnetic  field  intensity 

s ' 

H = far  scattered  magnetic  field  intensity 

r = distance  of  far  field  observer  from  origin 
o 

t/a  = normalized  time 
ka  = normalized  frequency 
a = characteristic  linear  dimension  of  target. 


The  scaling  and  normalization  that  is  indicated  above  yields  curves  which 
are  independent  of  target  size. 


The  impulse  response  augmentation  technique,  first  suggested  in 
1968  [2]  and  first  demonstrated  for  smooth  convex  targets  in  1973  [4],  and 
for  targets  with  edges  in  1974  [5] , deals  directly  with  the  smoothed  impulse 
response  of  the  targets  in  the  far  field.  The  smoothed  impulse  response  is 
computed  using  a space-time  integral  equation  approach  and  has  yielded  good 
results  up  to  body  sizes  of  several  pulse  widths  or,  equivalently,  up  to 
body  sizes  of  several  wavelengths.  The  regions  of  slow  variation  in  the 
smoothed  impulse  response  remain  the  same  in  the  exact  impulse  response; 
thus  it  is  only  necessary  to  determine  the  structure  of  the  singular  regions 
and  any  other  regions  of  fast  variation.  But  the  singular  portions  of  tlie 
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exact  impulse  response  that  result  from  scattering  by  specular  points  on 
smooth  convex  targets  can  be  computed  exactly,  an  i hence,  do  not  need  to  be 
computed  by  solving  the  space-time  integral  equation.  The  impulse  resr->onse 
augmentation  technique  combines  the  smoothed  impulse  response,  the  known 
singular  contribution  to  the  impulse  response,  and  the  theory  of  Fourier 
transforms  to  produce  the  total  impulse  response  and  the  frequency  response 
(system  function)  of  the  target  at  all  frequencies. 

The  impulse  response  augmentation  technique  is  most  easily  under- 
stood by  considering  the  most  basic  approach  to  the  deconvolution  (or  system 
identification)  problem.  A linear  system  (in  this  case  electromagnetic 
scattering  by  a target)  is  characterized  by  its  impulse  response  h(t)  or, 
equivalently,  its  system  function  (or  frequency  response)  H((jJ).  Of  course. 


h(t)  H(w) 

where  •->  denotes  Fourier  transform.  The  excitation  e(t)  of  the  linear  system 
in  this  case  is  the  regularized  (or  smoothed)  impulse 


e (t) 


n - (nt) ^ 
— e 

/rr 


(19) 


which  produces  the  regularized  (or  smoothed)  impulse  response  r(t)  of  the 
system.  This  response  is  given  by 

r(t)  = e(t)  * h(t)  (20) 


where  * represents  a convolution.  In  the  problem  being  considered  here, 
e(t)  is  specified  analytically  and  r(t)  is  computed  by  solving  the  space- 
time integral  equation.  It  is  desired  to  find  h(t)  and/or  H(gj).  This  is 
tlie  system  identification  or  deconvolution  problem. 


One  way  to  solve  tfiis  problem,  at  least  in  principle,  is  to  trans- 
form Eq.  (20)  and  rearrange  to  obtain 


H(U)) 


R(U)) 

E(w) 

f”'  1 11(0)} 


h (t) 
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(21) 


where 

E(OJ)  <-4  e(t) 

R(w)  ^ r (t) 

F Mh(w)}  is  the  inverse  Fourier  transform  of  H(w). 


However,  the  estimate  of  the  system  response  f(t)  that  is  computed  contains 
some  uncertainty  or  noise,  and  thus,  the  transform  of  the  computed  smoothed 
impulse  response  ft(CJ)  also  contains  some  noise.  In  using  Eq.  (21)  to  com- 
pute the  estimate  of  the  system  function,  it  can  be  shown  that  this  noise 
grows  exponentially  [4] , and  therefore,  this  brute  force  technique  will  not 
yield  valid  large  body  results. 

The  impulse  response  augmentation  technique  is  displayed  in  block 
diagram  form  in  Fig.  12.  This  technique  first  augments  the  smoothed  impulse 
response  to  remove  the  contribution  from  singular  portions  of  the  impulse 
response  that  are  known  exactly.  This  produces  the  augmented  smoothed  im- 
pulse response  r (t)  that  is  given  by 

r (t)  = r(t)  - e(t)  * f (t)  (22) 

a a 

where  f (t)  is  a suitable  augmentation  function  that  contains  the  known  sin- 

cL 

gular  portions  of  the  impulse  response. 


Next,  the  transform  of  the  augmented  smoothed  impulse  response, 

R ((jJ)  , is  computed  and  divided  by  the  transform  of  the  incident  pulse  to 

cl 

yield  the  augmented  frequency  response,  H^(w).  This  function  contains  noise 
which  increases  exponentially  at  frequencies  above  some  value.  However,  it 
is  known  that  the  augmented  frequency  response  must  go  to  zero  with  increa- 
sing frequency.  Thus,  an  estimate  of  the  high  frequency  behavior  of  the 
augmented  frequency  response,  ft  (w) , is  of  the  form 

3 


H (W)  ; W < 0) 

a c 


“a  " i 


F(C0)  ; h)  > U) 


(23) 
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FIG.  12  Impulse  response  augmentation  technique 


r 

I 

f 

I 


where  iy  t’ne  boundary  point  and  F(a))  is  the  hiyh  fr(‘f|uency  estimate  of 

(ah  . Tile  inverse  Fourier  tratisform  of  fi  (.;j)  tiien  yields  tiie  estimate  of 
o a 

of  the  augmented  irnijulse  response,  fi  (t)  . 

a 

Finally,  the  inverse  of  tlie  augmentation  jirocedure  is  performed  on 

/N 

h (t) , which  yields  the  estimate  of  the  impulse  res:x3nso,  h(t).  Moreover, 
a 

/V 

an  estimate  of  the  system  function,  H (w) , is  obtained  by  applying  the  inverse 

of  the  augmentation  procedure  in  the  frequency  domain  to  H (uj)  . 

a 

The  augmentation  function  represents  the  contributions  of  the  sin- 
gular portions  of  the  imfjulse  response  which  are  known  exactly.  These  sin- 
gular portions  may  contribute  to  not  only  the  high  frequency  behavior  but 
also  to  the  low  frequency  behavior  of  tlie  response,  as  in  the  case  of  an 
impulse.  These  contributions  are  removed  by  subtracting  the  effect  of  the 
augmentation  function  f (t)  from  the  response  to  yield  the  augmented  response 

cl 

as  given  in  Eq.  (22),  whicli  is  rei^eated  hero  for  convenience: 

r (t)  = r(t)  - e(t)  * f (t)  (22) 

a a 

or 

h (t)  = h(t)  - f (t) . 
a a 

Since  the  effect  of  the  high  frequency  contributions  has  been  removed,  then 
it  remains  to  estimate  the  manner  in  which  the  lower  frequency  components 
approach  zero  witii  incrcasinq  frequency. 

It  has  been  found  in  previous  work  [4-C]  that  the  augmentation 
functions  should  be  chosen  sucli  that  they  account  for  the  singular  contri- 
butions to  the  impulse  response  but  at  the  same  time  possess  a transform 
that  contains  only  linear  phase  variations.  Tlie  functions  wliich  satisfy 
this  simple  criterion  are  singularity  unctionals  and  pseudo-functions  that 
contain  only  a singl>'  dir.continuity . -Some  functionals  wiiich  possess  this 
characteristic  are  the  doublet,  the  impulse,  the  step,  and  Hadamard's  pseu- 
do-functions (10).  It  (la;;  also  been  found  for  smooth  convex  targets  that 
the  singular  portion  of  the  impulse  rcrg'onse  at  the  leading  edge  region 
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would  be  given  by  the  physical  optics  api>roximation  for  aspect  angles  where 
the  response  was  polarization  independent,  such  as  axial  incidence  on  rota- 
tionally  symmetric  targets.  For  the  case  of  smooth  convex  targets  with  a 
polarization  dependent  response,  the  first  order  physical  optics  correction 
will  yield  the  proper  singular  portion  of  the  impulse  response  at  the  leading 
edge.  For  the  case  of  targets  with  edges  the  singular  returns  from  the  edge 
regions  of  these  targets  is  obtained  by  use  of  the  results  obtained  from 
geometric  theory  of  diffraction. 

3.2  APPLICATION  TO  CYLINDER  WITH  FINS 

In  this  section  the  impulse  response  augmentation  technique  is 
applied  to  the  case  of  a cylinder  with  fins  attached.  This  work  extends  the 
techniques  that  were  developed  on  a previous  effort  [6]  for  the  case  of 
axial  incidence  to  the  case  of  oblique  incidence.  The  approach  used  in  this 
work  is  to  combine  the  two  tec)miques  previously  developed  for  open  thin 
surfaces  and  for  right  circular  cylinders  to  obtain  thie  impulse  response 
of  a cylinder  with  fins  attached. 

The  smoothed  impulse  response  of  this  target  was  obtained  using 
the  space-time  integral  equation  solution  technique  described  in  Sec.  2.2. 

The  target  is  centered  at  the  origin  with  the  axis  of  the  cylinder  coinci- 
dent with  the  z-axis  as  shown  in  Fig.  13.  The  cylinder  has  a length  to  dia- 
meter ratio  of  2:1.  Four  square  fins  with  a width  of  one  cylinder  diameter 
arc  placed  symmetrically  around  tlie  cylinder  body.  The  target  is  illuminated 
by  an  incident  plane  wave  witii  an  incident  pulse  width  equal  to  the  length 
of  tlie  cylinder.  The  smoothed  impulse  response  of  tliis  target  was  calcu- 
lated for  seven  angles  of  incidence  that  range  from  0°  to  180°  in  30°  incre- 
ments [5] . Tile  results  which  have  been  shown  to  be  in  good  agreement  with 
measurements  are  displayed  in  Fig.  1-1  for  0°  to  00°  and  in  Fig.  15  for  90° 
to  180°. 


The  impulse  response  augmentation  technique  next  re<}uires  a suit- 
able augmentation  function  as  discur.s.d  in  .^ec.  3.1.  For  the  cylinder  with 
fins  the  augmentation  function  will  be  taken  a;;  th<'  '.um  of  t 1r'  augmentation 
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incidence  for  TE  polarization. 
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P 


functions  for  the  rectunqular  plate  and  tlie  right  circular  cylinder  [6j . 
Thus,  the  aucjmentation  function  is  given  by 


f (t) 
a 


where 


D 

zi 


P 

2l 


S 

zi 


P 


z2 


M 


z3 


D 

Z 


D 5'(t/a-T  „ ) + P 6(t/a-T  ^ ) + S U(t/a-T  ) 
zi  zFi  zi  zFi  zi  zFi 


+ P 6(t/a  - T ^ ) + H (t/a  - T ^ - U(t/a  - T ) 

Z 2 ZF  2 Z 3 ZF  3 ZF  3 


+ D 6'(t/a-T  ) + P 6(t/a-T  ) 

z zC  1 z zc  1 


+ G (t/a  - T ) “3/2  ^ ) 

z ZC’  1 ZC'  1 


-1 

+ A (t/a  - T ) - U(t/a  - T ) 
1 ZC  1 ZC  1 


+ A (t/a  - T ) ‘ U(t/a  - T ) 
2 zc  2 zc  • 2 


t A (-t/a  t T ^ ) • U(-  t /a  + T ) 

3 zCi  ZC3 


+ f (t/a  - T ) 
a 4 zc  4 


the  cioublet  coefficient  given  by  t)ie  pliysical  optics 
approximation  for  the  fin  at  broadside  incidence 


the  impulse  coefficient  due  to  tlie  near  edge  of  ttie  fin 


the  step  coefficient  duo  to  the  response  from  tlie  sides 
of  the  fin 


the  impulse  coefficient  due  to  the  far  edge  of  the  fin 

-Ir 

the  t pseudo- function  coefficient  due  to  surface 
traveling  waves  on  the  fin 

the  (ioublot  coefficient  given  by  ttie  physical  optics 
approximation  for  the  cylinder  at  axial  incidence. 
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the  coefficient  for  the  cylinder  at  axial  inci- 

dence 


P 

z 


T 


zF2 


T 


zF3 


'^zCi''^zC2''’^zC3 


zC4 


f , (t/a  - T .) 
a4  zc4 


4 

B 


—3/2  . . 

the  t pseudo- f unct ion  coi'fficient  due  to  the  specular 

return  from  the  sid<?  of  the  cylinder  at  broadside  inci- 
dence 


the  t ‘ pseudo- function  coefficients  due  to  the  return 
from  the  edtji>s  of  the  cylinder 


the  time  at  which  the  leading  edge  of  the  impulse  resi>onse 
from  the  fin  begins 


tlie  time  at  which  the  singular  contribution  from  the  far 
eilge  of  the  fin  begins 


the  time  at  which  the  singular  contribution  from  the 
single  excursion  traveling  wave  on  tlie  fin  begins 


the  time  at  which  the  singular  returns  from  three  edges 
of  the  cylinder  begin 


the  time  at  whicli  the  cylinder  creeping  wave  appears 

-1  ^"'^'^zc4 

F A , e e 

4 


the  cylinder  creeping  wave  coefficient 


2.051 


The  values  used  for  the  above  parameters  are  summarized  in  Table  1. 
They  were  previously  obtained  for  the  case  of  the  right  circular  cylinder  and 
the  rectangular  plate  [5,6!. 


Once  the  augmentation  function  is  obtained,  the  augmented  frequen- 
cy response  is  calculated  by  subtracting  the  transform  of  tlio  augmentation 
function  from  tlic  frequency  rcisi'onse.  The  liigh  frequency  estimate  used  to 
obtain  the  estimate  of  the  total  augmented  frequency  response  takes  the  foim 
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TABLt:  1 


SUMMARY  OF  PARAMIOTEKS  USED  TO  Ri;PRr-;SENT 
THE  AUGMENTATION  FUNCTION  AND  THE  AUGMENTED  FREQUENCY  RESPONSE 


Angle  of  Incidence 


Parameter 

O 

0 

o 

c 

o 

60 

90" 

o 

120 

150° 

180° 

0 

0 

0 

1.273 

0 

0 

0 

0.701 

0.686 

0.955 

-0. 317 

0.  955 

0.686 

0.701 

s 

z 1 

-0.166 

-0.181 

-0.260 

0 

-0.260 

-0. 181 

-0.166 

‘’z2 

-0.224 

-0.049 

-0.318 

0 

-0. 318 

-0. 049 

-0.224 

“z3 

0 

0 

0 

0.203 

0 

0 

0 

D 

z 

0.5 

0 

0 

0 

0 

0 

0.5 

0.141 

0 

0 

0 

0 

0 

0.  141 

Ci 

Z 

0 

0 

0 

-0.318 

0 

0 

0 

-0.035 

0.  382 

0.290 

-0.123 

0.290 

0.382 

-0.035 

-0.284 

-0.016 

-0.081 

0 

-0.081 

-0.016 

-0.284 

0 

-0.106 

-0.012 

0 

-0.012 

-0. 106 

0 

'^zFl 

-4.0 

-3.464 

-2.0 

0 

0 

0 

0 

'^ZF2 

0 

0 

0 

0 

2.0 

3.464 

4.0 

n 

-a' 3 

0 

0 

0 

4.0 

0 

0 

0 

.T 

-4.0 

-4.464 

-3.732 
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where  the  cutoff  point  ui  is  chosen  sucli  tiiat  the  phase  of  li  (w)  is  still 

c a 


in  the  linear  reuion  and  A is  computed  so  that  H (a>  ) = H (W  ) . 

a c a c 


These  coef- 
ficients are  also  displayed  in  Table  1,  The  results  for  the  impulse  response 
and  the  frequency  response  are  displayed  in  Figs,  16-22  and  described  in  the 
following . 


3 . 3 RESULTS 


TE  impulse  response  of  the  finned  cylinder  for  axial  incidence  is 
displayed  in  Fig.  16a  and  the  corresponding  frequency  response  is  displayed 
in  Fig.  16b.  The  initial  return  of  the  impulse  response  at  t/a--4.0  con- 

-,V 

sists  of  a doublet,  impulse,  and  t " pseudo-function  from  the  cylinder  face 
and  edge.  Also  contained  in  this  leading  edge  of  tlie  impulse  response  is  an 

impulse  with  a positive  sign  due  to  the  fin  edge.  This  is  followed  at  t/a = 0 

— ^ 

by  a negative  impulse  from  tire  far  fin  edge.  At  t/a  = 4 there  appears  at' 
pseudo-function  from  the  far  cylinder  edge.  Finally,  a creeping  wave  can  be 
observed  at  t/a  ~ 6.  In  Fig.  16b  tire  corresponding  frequency  response  is 
displayed.  This  response  increa.ses  linearly  with  increasing  frequency  due 
to  tire  doublet  in  the  impulse  response  at  t/a  = -4.  Tire  ripples  in  this  re- 
-spoirse  are  mainly  due  to  the  interference  betweetr  this  doublet  and  the  re- 
turns from  the  far  fin  edges,  tire  far  end  of  the  cylinder,  and  the  creeping 
wave. 


The  TE  impulse  resfronse  and  frequency  r esponse  of  a finned,  cyl  iirder 

is  shown  for  a 30^  angle  of  iircidence  iit  Fig.  17.  The  initial  return  from 

the  near  edge  of  the  cylinder  face  appears  at  t '.i  -4.5.  This  is  tollowed 

at  t ""-3.5  by  a positive  impulse  aird  a nogativi'  step  from  tiro  iroar  fitr  edge. 

1 

Next,  a "backwarti"  t pr.eudo-f urrction  appears  at  t/.r~-2.5  from  the  f.rr 
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FIG.  16  TE  response  of  finned  cylinder  with  radius  a for  0°  angle 
of  incidence. 
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edge  of  tlie  front  face.  At  t/a = 0 a negative  impulse  is  observed  due  to  the 
far  edge  of  the  fin.  Following  this  is  a negative  t ""  pseudo- function  from 
the  far  edge  of  the  cylinder  at  t/a  =^2.5.  Finally  at  t/a -5.8  the  creeping 
wave  that  travels  around  the  cylinder  can  be  seen.  In  Fig.  17b  the  corres- 
ponding frequency  response  is  displayed.  This  response  approaches  a con- 
stant value  at  high  frequencies  due  to  the  two  impulses  from  the  edges  of  the 
fin.  The  slow  resonances  observed  at  the  higher  frequencies  are  due  to  the 
interference  between  these  two  impulses.  The  more  rapid  ripples  superimposed 
on  this  slow  waveform  are  due  to  interference  between  the  impulses,  the  t 
pseudo- functions  and  the  creeping  wave. 

In  Fig.  18  the  TE  impulse  response  and  the  frequency  response  of  a 
finned  cylinder  is  shown  for  a 60°  angle  of  incidence.  The  initial  return 
is  a t ^ pseudo-function  from  the  near  edge  of  the  front  face.  This  is  fol- 
lowed at  t/a  = -2  by  a positive  impulse  and  a negative  step  from  the  near  edge 

-A 

of  the  fin.  Next,  at  t/a  --.3  appears  a "backward"  t pseudo-function 

followed  at  t/a  = 0 bv  a negative  impulse  from  the  far  edge  of  the  fin.  Next 

_x 

at  t/a  -+.3  at"  function  occurs  due  to  the  near  edge  of  the  far  cylinder 
face.  Finally,  at  t/a  = 5.0  the  creeping  wave  tJiat  travels  around  the  cylin- 
der can  be  observed.  In  Fig.  18a  the  corresponding  frequency  response  is 
displayed.  As  the  frequency  increases,  this  response  continues  to  ripple 
about  a constant  value  due  to  the  interference  betvv^een  tlie  two  impulses  from 
the  fin  edges. 

The  TE  impulse  response  and  frequency  response  for  tiie  finned  cyl- 
inder at  a 90°  aspect  angle  is  sliown  in  Fig.  19.  The  initial  ret  irn  in  tlie 
impulse  response  is  a t pseudo-function  together  with  a negative  impulse. 

Next,  at  t/a = 0 appears  a doublet  and  a negative  impulse  duo  to  the  return 
from  the  face  and  edge  of  the  fin.  This  is  followed  at  t/a =3.1  by  the  ef- 
fect of  a creeping  wave  that  travels  around  the  cylinder.  In  Fig.  19b  the 
corresponding  frequency  response  is  shown.  One  can  observe  in  this  figure 
that  the  frequency  response  increases  linearly  with  frequency  due  to  the 
doublet.  Moreover,  it  oscillates  about  this  over  increasing  value  duo  to 
the  interference  between  the  doublet  and  the  t 


pseudo- func‘,  ion. 


HVkal  / 


Till-  Tt:  impulse  renpense  and  frequency  response  of  the  finned  cyl- 
inder is  shown  in  I'iq.  20  for  a 120  angle  of  incidence.  The  initial  return 
of  tile  impulse  response  is  a t ' pseudo -funct  ion  from  tiie  near  edge  of  tiie 
near  face  of  the  cylinder.  Noxt,  at  t/a--.3  a very  small  "backward"  t 

nseudo-Eunction  appears  from  the  Ear  edge  of  the  near  face  of  tiie  cylinder. 

-1 

An  impulse  due  to  the  near  fin  edge  is  observed  at  t/a = 0 followed  by  at" 
pseudo- funct ion  from  the  far  edge  of  the  cylinder  at  t/a  -+.3.  The  negative 
impulse  at  t/a  = 2 is  due  to  the  far  fin  edge.  Tiie  creeping  wave  at  t/a  = 5.0 
can  be  seen.  Tiie  corresi>onding  frequency  response  is  displayed  in  Fig.  20b. 
Again,  the  two  imj'ulses  ui  tiie  imiiulse  res[)onso  interfere  with  eacli  otlier  to 
give  an  oscillating  frequency  response  at  the  liigher  frequencies.  The  t ‘ 
pseudo-function  anii  tlie  creeping  wave  cause  distortion  of  what  would  otlior- 
v/ise  be  a very  regular  oscillation  in  the  frequency  resjjonse. 


Fig.  21  displays  the  TE  impulse  response  and  frequency  response  of 

a finned  eviinder  for  a 15o”  angle  of  inci'iencc.  The  initial  part  of  this 

_ 1 

response  at  t/a  = -4 . 5 is  a t " pseudo-function  from  the  near  edge  of  the  near 
face  of  the  cylinder.  This  is  followed  at  t/a =-2.5  by  a "backward"  t " 
pseudo- funct ion  from  the  far  edge  of  the  near  face  of  the  cylinder.  At 
t/a = 0 the  impulse  and  negative  step  from  the  near  fin  edge  occui's.  Follow- 
ing  this  at  t/a  =2.5  at"  pseudo-function  appears  due  to  tlie  far  edge  of 
the  far  cylinder  face.  This  is  followed  b'/  a small  negative  impulse  at 
t/a = +4  from  the  far  fin  edge.  Finally,  the  creeping  wave  can  be  seen  at 
t/a  =5.8  in  the  impulse  response.  The  corresponding  frequency  response  shown 
in  Fig.  21b  again  approaches  a constant  value  at  high  frequencies.  Tiie  rip- 
pling in  ttiis  structure  is  due  to  the  interference  between  tiie  impulse,  tiie 
- i- 

t " pseudo- funct ions,  and  the  creeping  wave. 


Fig.  22  ilisplays  the  TE  impulse  response  and  freguenev  rostjonse 
of  a finned  cylinder  for  a 18o”  an<jle  of  incidence.  Tiie  doublet  and  impulse 
that  appear  at  tlie  leading  edge  of  tiiis  response  are  duo  to  tiie  si  ecular  and 
edge  returns  from  tlie  front  fact'  of  tlie  cylinder.  l.ie  imi'ulso  and  negative 
step  at  t/a  = 0 are  due  to  the  return  from  tlie  near  fin  edge.  Tlie  return 
duo  to  the  far  fin  edge  and  far  cy'linder  ai'pears  at  t/a  = •!  and  finalK'  tlie 
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FIG.  21  TE  response  of  finned  cylinder  with  radius  a for  150  angle 
of  incidence. 


creeping  wave  can  be  seen  at  t/a=^6.  The  corresponding  frequency  response  is 
displayed  in  Fig.  22b.  This  frequency  response  increases  linearly  with  in- 
creasing frequency.  Furthermore,  it  oscillates  about  this  linearly  increasing 

-A 

value  due  to  the  interference  between  the  two  impulses,  the  doublet,  the  t ■' 
pseudo- functions , and  the  creeping  wave. 


-54- 


SKCTXON  4 


TTMK  DOMAIN  INVERSE  SCATTEKINti  - AXIAL  INCIDENCE 

Tho  invorso  acatteriaq  problem  consists  of  determining  tlie  struc- 
ture of  .111  unknown  scutterer  given  information  about  a field  incident  upon 
the  scutterer  and  the  resi)onse  of  tlic  scatterur  to  that  field.  Most  of  the 
uttempis  at  solving  the  inverse  scattering  problem  have  been  carried  out  in 
tlie  treguency  domain.  The  inherent  advantages  of  employing  time  domain  tech- 
nitiuec  appear  to  have  been  first  docuinenteti  by  Kennauqii  and  Cosgriff  [Hi. 
Tiit'Y  show<.'d  that  if  phy.sical  o|)tics  currents  were  postulated  on  the  .surface 
of  the  scutterer,  then  tliey  jiroduced  an  approximate  impmlse  resiJonso  that  was 
simply  tlie  seconti  derivative  of  the  projected  area  function  of  the  scatterer. 

Many  researchers  liave  devcloxied  approaches  [12,13i  to  the  inverse 
scattering  problem  tiiat  liavc  been  based  on  extensions  of  tliis  result.  An- 
othe''  approach  based  on  polarization  information  in  tlie  imj->ul.se  ro.si  onse  lias 
also  been  dcvelopei-i  I14,.15].  These  approaches  have  reiiroduccd  the  scatterer 
geometry  with  varying  degrees  of  success.  This  would  be  expected  since  the 
relation  between  the  impulse  response  and  the  two  derivatives  of  the  area 
function  is  exact  only  at  the  leading  edge  of  the  scattered  field  response, 
a single  point  in  time.  After  the  leading  edge,  the  response  is  altered  by 
currents  arriving  from  other  space  points.  Tlicrcfore,  the  ;)hysical  optics 
solution  must  be  "corrected"  by  tiiese  currents  flowing  on  the  body.  Eor  a 
given  object,  if  tlie  incident  pulse  width  is  short  compared  tc>  body  size 
(the  optical  limit)  then  the  correction  currents  will  have  a small  eif'-.  t, 
aiui  optical  ray.s  can  tit'  |ilaced  in  oiu'-to-one  correspon  ience  with  point  s on 
til’-'  body.  On  the  other  liaiui,  if  the  size  of  the  Ixi  ly  is  comparabK'  to  a 
puls.c  wi'itli,  then  thi>  "correction"  terms  have  a stiong  effect  on  the  ■inti  ':: 
and  the  physical  of'tics  solution  Is  di"U'aded.  In  the  case  s-m.i  1 1 , 

tlie  correction  terms  liominate  the  result  and  the  pliy.sical  opto  --,  s.dntion  i - 
noaningles 

In  this  stu  iy,  the  invetse  scattering  i rolden  is  tormulat'-.i  ,i.  .u, 
invers.ion  of  the  sp.ice-t  ime  integral  cguation.  Ttiis  ai's-oach  \o  ■ 1 ■ i n i n iiii 
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j 

i 

I 
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the  response  of  the  scatterer  incorporates  both  the  effects  of  the  incident 
field  directly  and  those  due  to  the  correction  currents  flowing  on  the  body 
and  was  first  suggested  in  1974  [5] . By  using  a complete  description  of  the 
surface  current  interactions,  this  formulation  should  yield  a very  close 
approximation  to  the  target  geometry. 

4.1  DERIVATION  OF  EQUATION.S 

In  order  to  derive  the  inversion  procedure,  it  is  necessary  to 
first  examine  the  direct  solution  of  the  scattering  problem.  The  direct 
problem  consists  of  determining  the  scattered  field  given  a description  of 
the  incident  field  and  the  scatterer  geometry.  This  problem  is  solved  by 
first  computing  the  currents  flowing  on  the  scatterer  surface.  From  these 
currents  the  scattered  field  can  be  calculated  directly. 

The  expression  for  the  surface  currents  is  derived  from  an  expres- 
sion for  the  total  field  at  an  arbitrary  point  in  space  (see  Sec.  2.2.1) . 
This  arbitrary  point  is  specialized  to  a point  on  the  scatterer  surface  and 
the  appropriate  boundary  conditions  are  applied.  This  yields  an  equation 
for  the  surface  current  J at  point  r on  the  scatterer  surface  and  time  t as 

J(r,t)  = 2a  X H^(r,t)  + J (r,t)  (24) 

n c 

where 


J (r,t)  = 
c 


_1 

2tt 


H (r,t)  = incident  magnetic  field 


a = unit  vector  normal  to  surface 
n 


r = position  vector  to  observation  [Joint 


r'  = position  vector  to  integration  point 
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K = i - r 


r - r 


t Lime  in  1 1 qiit -m. ■ t' -r : 


The  firut  term  on  the  riqliL-liami  nuie  at  hq.  (24)  may  be  considered 
thc>  source  term  and  represents  tlie  direct  influence  of  the  incident  tield  on 
the  current  at  the  observation  point  (r,t).  Moreover,  tliis  term,  wlien  ap- 
plied to  the  illuminated  side  of  the  scatteror,  yields  the  familiar  physical 
optics  approximation  for  the  surface  current.  The  second  term  on  the  right- 
hand  side  of  Eg.  (24)  represents  the  influence  of  currents  at  other  surface 
points  on  the  current  at  (r,t).  It  is  important  to  note  that  the  influence 
of  the  current  at  other  points  on  the  surface  on  the  current  at  (r,t)  is 
delayed  in  time  by  R,  the  distance  between  the  two  ixiints.  This  allow's  Eg. 
(24)  to  be  solved  by  a "marching  on  in  time"  procedure  rather  than  necessi- 
tating matrix  inversion. 

Once  the  surface  currents  liave  been  determined  the  far  scattered 
field  can  be  calculated  from  them  using 


->-s  ->■  1 

H (r,t)  = 


4iir  9t 
o 


J ( r ' , 1 ) X a 


dS 


(25) 


T = t - K 


wliere 


r = distance  to  the  far-fiold  observer 
o 

a = unit  vector  from  tlic  integration  point  to 
^ the  far-fiold  observt'r. 

If  the  surface  current  oxpres.sion  in  Eg.  (24)  is  substitutcci  in  Eg.  (25), 
tlien  tile  result  is 
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~^s  1 3 ™ ^ ^ 

V‘  .1  2a^  X H (r  ,T)  X a^  ds' 


T = t - R 

■*■7^^  ^ (r",T)  X a ds'  (26) 

4~  oT  .1  c r 

T = t - R 

If  the  incident  field  is  an  impulse,  then  the  first  term  in  Eq.  (26)  is  sim- 
ply the  term  that  has  been  recognized  previously  [11]  to  be  proportional  to 
the  second  derivative  of  the  area  function.  Assuming  the  incident  wave  is  an 
electromagnetic  impulse,  then  Eq.  (26)  becomes 


- 1 1 3 - ^ « 

" 2i  -^7-  “h  * « 57  J ■’cl*'  '">  " *r 


where 


T = t - R 


->S 

r^H^(r,t)  = the  impulse  response  of  the  target 


S(t  ) 
s 


the  silhouette  area  of  the  scatterer  as 
delineated  by  the  incident  impulse  assumed 
moving  over  the  scatterer  at  one  half  the 
free-space  velocity 


r = distance  of  far-field  observer  from  the  origin 
o 


t = t + r 
s o 


” Ih'i 


J = J that  results  from  an  incident  impulse, 
cl  c 


This  equation  may  be  simplified  by  integrating  twice  to  obtain 
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(28) 


wiiero 


■'■s  ->•  1 . 1 

r li  (r,t)  = — - F;(t  ) a + — ttt  ' J (r  ,T)  x a 

ok  2 II  s 'it  . cK  r 


T = t - k 


>s 

= the  ramp  response  of  tlie  target 


J = J tiiat  results  from  an  incident  ramp  waveform. 
cK  c 


Thus,  by  direct  consideration  of  the  space-time  integral  equation  the  exact 

relationship  between  the  target  response  and  the  target  geometry  has  been 

obtained.  In  particular,  Eq.  (28)  gives  the  target  ramp  response  in  terms 

of  both  the  target  area  function  and  the  contribution  due  to  t.he  "correction 

currents,"  J . Moreover,  it  is  important  to  note  that  the  corr._-ction  cur- 
c 

rents  as  given  in  Eq.  (24)  are  time-retarded  functions  of  currents  at  other 
space  points,  and  thus  will  be  zero  at  the  leading  edge  of  t)ie  incident  wave- 
front  as  it  travels  across  the  target.  It  is  this  feature,  exclusive  to  tlie 
time  domain  formulation,  that  allows  tlie  determination  of  tlie  target  geometry 
from  its  ramp  (or  equivalently,  impulse)  response. 


4.2  NUMERICAL  SOLUTION 

In  order  to  test  the  effectiveness  of  this  approach,  the  technique 
was  applied  to  the  class  of  rotationally  symmetric  scattorers.  T)ie  rota- 
tionally  symmetric  scattering  problem  is  depicted  in  Fig.  23.  For  this  case 
the  scatteror  is  symmetric  about  the  z-axis,  the  incident  field  is  axially 
incident,  and  the  far  field  is  computed  in  the  baclcscatter  direction.  Tlie 
contour  of  a rotationally  symmetric  object  can  be  completely  described  by 
tlie  radius  vector  p which  varies  as  a function  of  z,  and  its  ;rojected  area 
function  can  be  expressed  simplv'  as 

S = TTp^  (z,t)  . 
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Subb’tituting  this  projected  area  function  into  Eq.  (28)  and  rcar- 
ranqinq  terms  yields  the  inversion  equation  for  the  rotationally  symmetric 
case 


P(z,t)  = 


2r  H^(r,t) 
o R 


JL  _il 

2TT  3t 


J ( r ' , T ) X a 
cR  r 


ds' 


(29) 


T = t 


This  equation  gives  p(z,t)  in  terms  of  tiio  target  ramp  response,  which  is 
known,  and  in  terms  of  the  correction  currents  at  earlier  times,  wiiich  have 
eitJier  been  previously  computed  or  are  known  to  be  zero. 

To  solve  this  equation  a technique  was  developed  under  a previous 
program  [5]  which  iterated  on  estimates  of  the  entire  target  geometry.  For 
this  program  it  was  proposed  to  improve  this  technique  using  a direct  time 
domain  approach.  In  this  approach  the  target  geometry  is  generated  sequen- 
tially in  time  as  the  incident  field  moves  across  the  target.  Geometry 
values  are  determined  from  the  incident  field,  previously  computed  geometry 
values,  and  correction  currents  set  up  on  the  structure  as  a result  of  sur- 
face current  interaction. 


A flow  chart  of  the  direct  solution  procedure  is  shown  in  Fig.  24. 
The  general  ai'proach  is  to  step  along  the  target  contour  in  equally-spaced 
intervals.  Tlie  value  of  p(z,t)  is  computed  at  the  center  of  each  interval 
at  the  time  at  whicfi  the  initial  return  from  that  point  readies  the  far 
field  observer.  This  insures  that  the  correction  current.s  at  that  time  are 
a function  only  of  currents  flowing  on  the  portion  of  the  target  who.se  con- 
tour is  already  know-n. 

The  sol  .lion  procedure  beciins  witli  an  estimate  of  . (z,t)  for  tiu 
first  interval  along  the  target  contour.  The  jdarcmeiit  of  tlie  inte  rval  i.. 
obtained  by  assuming  the  target  is  locally  spherical  at  t!ic  specular  [xiint. 
The  value  of  p(z,t)  at  the  center  of  tlie  interval  is  computed  from  Ei.  (J9) 
by  neglecting  the  surface  integral,  giving 

*•<;  ->■  i 

0 (z  , t ) - 1 2r  II’  (r  , t)  1 " . 

1 o R 
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From  this  estimato,  the  qcomotry  parameters  necessary  for  the  calculation  of 
the  correction  currents  are  structured.  This  is  accom^^dished  by  constructing 
a band  around  the  target  whose  width  is  equal  to  that  of  the  contour  inter- 
val. The  band  is  then  divided  into  patches  of  equal  area.  Tlic  location  of 
tlie  center  cf  tlie  band  z.,  tlie  corresponding  radius  vector  p(z^,t),  and  tlie 
components  of  the  unit  normal  at  p(z  ,t)  are  also  necessary  inputs. 

i 

Once  tiie  geometry  of  the  first  band  has  been  structured,  tiien  an 
estimate  is  made  of  the  location  of  tiie  center  of  the  next  band.  This  is 
obtained  by  projecting  the  unit  normal  at  p(z^,t)  forward  from  p(z^,t)  for 
one  interval  length  and  then  finding  the  z^  corresponding  to  the  end  of  the 

interval.  In  order  to  compute  the  value  of  p at  z^/  it  is  necessary  to  com- 
pute the  correction  currents  at  tlie  time  t^  at  which  the  incident  field 

reaches  z^.  Tliis  computation  is  done  by  an  improved  version  of  the  program 
ROTSY  which  solves  the  direct  scattering  problem  for  rotationally  symmetric 
targets.  The  details  of  this  computation  are  described  in  Appendix  8.1. 

The  radius  vector  p(z^,t)  is  then  computed  from  Eq.  (29)  using  the  correction 
currents  at  t^  and  the  value  of  t.he  far  scattered  field  -Jt  the  time  t^^ 
which  the  initial  return  from  z^  reaches  the  far  field  observer. 

The  accuracy  of  tlie  results  obtained  from  the  ROTSY  program  is 
improved  if  the  target  is  divided  into  approximately  equal  bands.  In  order 
to  insure  this  for  the  inversion  solution,  a second  estimate  is  made  of  the 
location  of  the  band  center,  2.^.  This  estimate  is  obtained  by  constructing 
the  band  along  the  line  between  P(z^,t)  and  the  first  estimate  of  n(z^,t). 
.Again,  the  correction  currents  are  computed  for  the  second  estimate  of  z^ 
and  p(z^,t)  is  recalculated  from  Eq.  (29). 

Each  value  of  p(z.,t)  is  tested  to  determine  if  it  i .s  negative. 

1 

If  it  is,  then  the  inversion  procedure  terninatos.  If  it  is  not,  then  the 
geometry  {laramcters  for  band  1 are  computed  and  the  structuring  of  band 
i t-  1 begins.  In  this  manner,  the  inversion  technique  generate.s  the  target 
geometry  sequentially  in  time. 
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4 . 3 RESULTS 


In  order  to  test  t!ie  direct  solution  procedure  described  in  the 
previous  section,  a computer  proyram  INVERD  was  written  to  implement  the 
{-■'rocedure . This  program  was  used  to  compute  the  contours  of  the  four  rota- 
tionally  symmetric  objects  shown  in  Fig.  25: 

(1)  A spltere 

(2)  A cylinder  witli  two  sphere  caps 

(3)  A cylinder  with  a flat  front  end  and  a sphere-capped 
back  end 

(4)  A cylinder  with  a si>here-capped  front  end  and  a flat 
back  end. 

For  the  four  objects,  the  radius  of  the  spherical  segments  was  0.5  meters. 

The  radius  of  the  cylinder  body  was  also  0.5  meters  and  its  length  was  1 
meter.  In  all  cases  the  incident  field  was  axially  incident  and  the  far 
field  was  calculated  in  the  backscatter  direction.  The  ramp  responses  of 
the  four  objects  are  depicted  in  Fig.  26.  Since  these  same  four  objects 
wore  also  used  to  test  the  iterative  solution  technique  which  was  previously 
developed  [5] , valid  compiarisons  can  be  made  between  the  two  approaches. 

The  implementation  of  the  direct  solution  procedure  gave,  in  gen- 
eral, very  good  results.  For  all  four  targets  the  contours  were  reproduced 
with  a high  level  of  accuracy.  The  only  discrepancies  tended  to  occur  at  t)ie 
back  of  the  targets  whore  the  target  size  was  slightly  overestimated.  Per- 
spective plots  were  made  of  the  actual  target  contours  and  the  ap£>roximations 
to  the  contours  using  both  the  direct  and  iterative  procedures.  Tlic  coordi- 
nate axis  and  the  view  angle  for  the  plots  is  shown  in  Fig.  27. 

The  results  obtained  for  the  four  targets  are  shown  in  Figs.  28  - 
31.  For  the  case  of  the  sptiero  (Fig.  28),  tlic  results  are  very  close  to  the 
actual  contour.  Tlie  front  of  the  spliero  is  rej 'reduced  exactly  and  there  is 
only  a slight  overest imaf ion  of  the  back.  For  the  cylinder  with  two  s)'here 
caps  (Fig.  29),  again  ttrero  is  the  exact  reprociuct  ion  of  Uu:  front  ;:phore 
cap  and  the  sliq)it)y  enlarged  back  sph-'re  ca;  . The  cylinder  side  is  enti- 
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(ii)  Iterative  Solution 


(bl  Direct  Solution 


(cl  Actual  Contoui 


FIG.  28  Contour  estmiiitos  for  n spherr?. 
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1 


I 


(c>  Actual  Contour 


FIG,  29  Contour  nstimates  for  a cylinder  with  two  sphere  caps. 
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(a)  Iterative  Solution 


(c)  Actual  Contour 


FIG.  30  Contour  estimates  for  a cylinder  with  a flat  front  end  and 
a sphere  capped  hack  end. 
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(c)  Actual  Contour 


FIG.  31  Contour  estimates  for  a cylinder  with  a sphere- capped 
front  end  and  a flat  back  end. 
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mated  quite  accurately  and  is  also  somewhat  smoot)ier  than  estimate  obtained 
using  the  iterative  procedure. 

Fig.  30  shows  the  results  for  the  cylinder  with  the  flat  front  end 
and  the  sphere-capped  back  end.  This  estimate  also  very  closely  approximates 
the  target  geometry.  The  flat  front  end  is  accurately  recovered  and  there  is 
again  less  oscillation  along  the  cylinder  sides  than  for  the  iterative  re- 
sults. As  with  the  previous  targets,  there  is  the  overestimation  of  the  back 
sphere  cap.  It  should  be  pointed  out  that  for  this  target,  it  was  deter- 
mined from  the  initial  p(z,t)  estimate  that  the  front  end  of  the  target  was 
flat  and  the  geometry  parameters  were  structured  accordingly. 

The  results  for  the  cylinder  with  the  sphere-cap  front  and  flat 
end  back  (Fig.  31)  were  not  as  accurate  as  for  the  other  targets.  This  was 
expected,  however,  because  of  the  difficulty  in  obtaining  the  numerical 
accuracy  in  the  correction  current  computation  necessary  to  exactly  cancel 
the  tail  of  the  ramp  response.  Nonetheless,  the  direct  approach  did  give 
a fairly  close  approximation  to  the  flat  back.  It  also  very  accurately  re- 
produced the  sphere-capped  front  and  cylinder  sides. 

In  comparing  the  direct  solution  to  the  iterative  solution,  the 
results  indicate  that  the  direct  approach  gives  a better  estimate  for  some 
cases  while  the  iterative  approach  is  better  for  others.  The  iterative  solu- 
tion is  more  accurate  in  reproducing  the  back  sphere  cap  while  the  direct  so- 
lution tends  to  approximate  the  cylinder  sides  more  effectively  and  in  gen- 
eral, give  smoother  target  contours.  The  direct  solution  does,  however, 
have  one  very  important  advantage.  The  solution  teclwique  allows  for  a sub- 
stantial decrease  in  run  time.  The  run  time  using  the  direct  approach  is 
cut  by  a factor  equivalent  to  the  number  of  iterations  necessary  for  con- 
vergnneo  in  the  iterative  technique.  This  advantage,  combined  with  the 
rather  marginal  differences  in  the  results  obtained,  makes  tlie  direct  tech- 
nique a more  favorable  approacli. 
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SECTION  5 


f 


TIME  DOMAIN  INVERSE  SCATTERING  - OBLIQUE  INCIDENCE 


In  thin  section  tlic  problem  of  inverse  scattcrin'j  by  rotational!'/ 
symmetric  targets  witii  obli(|ne  incitltmce  is  considered.  It  is  siiown  that 
tlie  geometry  of  a scatterer  can  be  obtained  from  simultaneous  processing  of 
its  responses  at  t'.vo  po lar  i za t i on.s  for  the  class  of  rotationally  syinmotric 
targets  at  arbitrary  (known)  aii'jle  of  observation.  .Scctie.n  5.1  gives  a 
baclcqround  of  the  problem  and  the  inverse  scattering  oijuations  are  derived 
in  .Sec.  5.2.  The  first  order  rami'  response  of  a curvilinear  square  patch 
is  developed  in  Sec.  5.3  and  used  in  the  numerical  solution  procedure  that 
is  described  in  Sec.  5.4.  The  results  for  three  targets  are  presented  in 
Sec . 5.5. 

5 . 1 GENERAL  PROBLEM 

It  has  been  known  for  some  time  that  the  impulse  response  of  a 
scatterer  in  the  pliysical  optics  limit  is  the  second  derivative  of  its  pro- 
jected area  function  [11].  Equivalently,  tlie  ramp  response  itself  is  simply 
the  projected  area  function.  Based  on  this  ph'/sical  optics  appro.ximation , 
reasonably  accurate  schemes  have  been  developed  for  obtaining  target  geom- 
etry from  the  area  function  [12,13].  This  technique  is  strictly  valid  only 
in  the  limit  of  vet','  high  frequencies,  or  equivalently,  at  t'ne  leading  edge 
of  the  rami''  response.  When  the  wavelength  is  of  the  order  of  tlie  target 
dimension,  the  solution  is  degraded  due  to  the  existence  of  surface  current 
interactions.  Using  the  space-time  integral  equation  tociiniques,  the  ef- 
fects of  these  "correction  currents"  can  be  included.  This  was  fir.st  done 

[5]  for  rotationally  symmetric  bodies  at  axial  incidence,  usinu  an 
iterative  teciiniquc.  In  Sec.  4 of  tliis  repiort  th'  technique  was  moiiified 
for  direct  solution  for  this  class  of  targets. 
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n tiuj  piiysical  optics  approximation,  thu  response  is  polariza- 
tion 1 ndt'i'undent . Wlien  t’ne  effect  of  correction  currents  is  included  in 
the  aiialys  is , it  is  seen  tiiat  tlie  difference  in  the  resjionses  at  polariza- 
tions in  two  directiotis  is  related  to  t!ie  difference  in  tlie  surface  curva- 
tures in  two  directions.  Tliis  fact  permits  treatment  of  a more  general 
class  of  targets  provided  tiiat  tlie  relationship  between  local  curvature  and 
overall  surface  contour  is  known.  In  this  report  it  is  assumed  that  the 
surface  is  convex  and  rotationally  symmetric  ana  tiiat  tiie  angle  of  incidence 
witii  respect  to  the  axis  of  rotation  is  known.  In  addition  symmetry  about  a 
plane  perpendicular  to  tlie  axis  of  rotational  symmetry  is  assumed.  It  is 
believed  tiiat  tlie  removal  of  tliese  (or  equivalent)  symmetry  conditions  from 
the  general  inverse  scattering  problem  requires  observation  at  several  di- 
rections. A sketcii  of  tills  target  class  appears  in  Fig.  32. 

While  the  tliooretical  development  which  follows  is  valid  for  tlie 
class  of  targets  described  above,  tlie  numerical  implementation  of  tlie  tech- 
nique was  limited  to  broadside  incidence  (i.e.,  perpendicular  to  the  axis 
of  symmetry).  The  tost  cases  were  a sphere,  a prolate  spheroid  and  a cylin- 
der. 

5.2  DERIVATION  OF  INVERSE  SCATTERING  EQUATIONS 

We  write  the  far-fiold  response  in  terms  of  the  time-retarded  sur- 
face currents: 


where 


r H(r  , 
o o 


"f> 


Jl  _1 
4t:  dt 


f 

s 


J(r,t) 


X a dS , 
r 

o 


(30) 


a IS  the  direction  of  observation, 
r 

o 


t = t-  + a • r , 
f r 

o 


r is  a point  or,  the  surface  S. 
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Tl'io  time  oriqins  aru  chosen  such  that  an  impulse  from  the  origin  at  t = 0 
arrives  in  the  far  field  at  t^ = 0.  The  units  of  time  are  light-meters.  We 
will  divide  the  surface  current  at  a point  into  two  parts,  that  due  to  tiie 
incident  field  and  the  part  duo  to  the  otiier  surface  currents: 


J(r,t)  = J (r,t)  t J (r,t) 
Po  c 


(31) 


\;'nere 


J (r,t)  = 2a  X H (r,t)  (32) 

Po  n 

the  part  due  to  tlie  incident  field,  is  identical  to  the  surface  current  in 
the  zero  order  physical  optics  approximation;  and 


J(r,t)=7^  j a x[_j-J(r',T)xaldS  (33) 

c 27!  J n R 

S 

is  called  tlie  "correction  curreiit".  liero, 

1 1 3 

f IS  the  operator  — + — , 

r''  ^ 

r”  is  a point  on  the  surface. 


a,^  = (r-r')/R, 

= t - K , 

a i::  the  :-;ut  face  normal  at  r . 
n 

The  incident  field  is  assumoi  to  b<-  a r.imi  arriving  at  t)i>'  origin  at  t ' 0: 


i 
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~^inc  ^ ^ ^ -*■ 

H (r,t)  = a„(t  + a * r)  , for  (t  + a • r)  >0 

Hr  r 

o o 


(34) 


= 0 


, for  (t  + a • r)  < 0 
^0 


The  solution  of  these  equations  will  proceed  in  the  usual  method  of  step- 
ping along  in  time.  Accordingly  we  divide  S (the  surface  illuminated  at  the 
retarded  time  t corresponding  to  far  field  time  t^)  into  a part  which  is 
presumed  )cnown  from  determination  at  an  earlier  instance  of  t^  and  the  new 
part  As,  which  is  unknown  and  which  is  just  illuminated  by  the  beginning  of 
the  incident  ramp.  Hence,  we  rewrite  Eg.  (30)  as 


r H 
o 


'"o'"f>  = 


_1  _3 
41T  3t 


J (r,t)  X a 


As 


dS  + 


4tt  3t 


J (r,t)  X a dS  (35) 


S-As 


Since  the  incident  ramp  has  just  passed  over  the  new  area  As  there  are  no 
contributions  from  AS  to  the  surface  currents  J(r,t)  for  r in  (S-As) ; i.e., 
the  integral  in  Eq.  (33)  is  also  over  the  known  area  (S-As) . Similarly,  at 
a point  r in  AS,  there  can  be  no  contribution  to  J(r,t)  due  to  surface  cur- 
rent in  (S-As) . The  two  integrals  in  Eq.  (35)  are  thus  completely  separa- 
ted. 


Consider  first  the  integral  over  the  known  region.  To  evaluate 
this,  the  surface  current  J must  be  computed.  The  integral  Eq.  (33)  con- 
tains  r as  an  integration  point.  It  is  necessary  to  separate  out  the  effect 
of  the  local  region  S^  about  r in  order  to  obtain  J(r,t)  explicitly  in 
terms  of  the  incident  field  and  the  retarded  currents  elsewhere  on  the  sur- 
face. Analytically  S can  be  considered  to  be  an  infinitesimal  patch  about 
>■  -* 
r,  while  in  the  numerical  solution  S^  is  the  surface  patch  with  center  r 

of  the  numerical  representation  of  the  target.  The  integration  over  is 
performed  in  Appendix  8.3  and  gives  rise  to  a "self  {latcli  correction  fac- 
tor", Y which  is  dependent  upon  polarization,  patcli  curvature  atid  patcli 
size.  In  accordance  with  tlie  result  of  Appendix  8.3  we  reformulate  ex- 
pressions (31) -(33)  as 
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(36) 


where 


'j(r,t)  = J (r,t)  + J (r,t) 
POl  c 


is  the  first  order  physical  optics  current  and  includes  the  effect  of  the 
self-patch;  and  where 


->■->•  1 f» 

J (r , t)  = — — 
c 2TrY  j 


a X [(■  J (r  ^ , T)  X a ] dS 
n 


is  the  correction  current  due  to  the  remainder  of  the  surface.  Note  that 
J in  this  integral  in  turn  consists  of  the  sum  of  the  retarded  first  order 
physical  optics  currents  and  the  correction  currents  at  the  integration  point 

-V  ^ -> 

r , as  in  Eg.  (35) . In  the  numerical  solution  we  need  only  store  J at  all 

points  (r  ,T)  and  compute  the  contribution  due  to  directly  when 

integrating  Eq.  (38).  Also,  as  will  be  indicated  later,  the  quantity  of 

interest  is  actually  9J  (r,t)/9t.  Differentiating  Eq.  (38)  and  integrating 

c 

numerically  over  the  discrete  patches,  we  can  take  the  derivative  inside  the 
integral  for  these  patches  which  are  fully  covered  by  the  incident  field, 
while  we  need  special  treatment  for  the  patches  which  are  only  partially 
excited  by  the  incident  field  and  for  which  9s/9t  is  non-zero.  Anticipating 
the  results  of  the  next  subsection  we  will  approximate  tliis  effect  by  sub- 
stituting an  effective  area  for  the  nominal  patch  area. 

For  an  incident  ramp  we  have  simply: 

9 ;^inc,“>’'  1 

(r  ,T)  = --  a^^ 

R 

for  (T  + a • r ) >0  and  zero  otherwise.  We  thus  write  numerically: 
r 

o 


9j  (r,t) 
c 


9t  2tiy  . n 


ax  V 


^'Villi’ 

--  + J ^ 


non-selt  » x. 
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Where  the  summation  is  over  all  patches  £ except  the  self-patch  and  where 

is  understood  to  be  an  effective  area  for  partially  excited  patches. 

This  is  the  expression  that  is  coded  in  the  computer  program.  The  process  of 

evaluating  this  expression  will  be  referred  to  later  as  the  process  of  "up- 

dating  the  correction  currents".  Note  that  J at  time  t is  given  in  terms 
->  ^ 
of  the  J at  other  points  at  previously  calculated  time  point  T = t - R. 
c 

We  are  now  in  a position  to  compute  the  second  integral  of  Eg.  (35) 
to  obtain  the  contribution  of  the  >cnown  area  (S  - AS)  to  the  far  field.  As 
was  done  above,  we  take  3/9t  inside  the  integral,  keeping  in  mind  the  diffi- 
culty of  partially  excited  patches.  Substituting  Eq.  (36)  into  the  second 
integral  of  Eq.  (35)  we  write: 


II 


(S-AS)  " VpOI^^o'^^  ^ ^o“c‘^o'^f^  ' 


where 


rH_^(r,t)=—  ) — a.x-r—  J (r,t) ; x a As,  , 

o POl  of  4tt  Y.  , ni  9t  PO  , 

3 3' 


"o  ^ 


or,  making  use  of  Eq.  (5)  and  expanding  the  vector  triple  product: 


"o%l^^o' 


t^ ) = -^  Y — a . *a  ) As.  , 

f 2T7  / . Y . n]  r 3 

3 3 o 


where  summation  is  over  all  patches  with  r^  in  (S-AS) . We  recognize  this  as 
the  projected  area  function  of  the  physical  optics  approximation,  modified 
by  the  individual  self-patch  correction  factors,  Y ^ . (We  also  remember  that 
As.  is  an  effective  area  for  partially  excited  patches.)  The  other  part  of 
II  is  written  as 


9j 

r H (r  ,t  ) = — (r,t)  x a dS 

o c o f 4tt  ot.  r 

. o 

S-As 


(40) 


This  is  also  evaluated  numerically  by  summing  over  all  patches  with  r^  in 
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r 


(S-AS)  . (In  this  integral  we  need  not  '.'orry  about  partially  excited  patches, 
since  the  correction  currents  are  still  zero  at  that  time,  so  that  such  pat- 
ches can  simply  be  omitted.) 


We  return  now  to  the  first  integral  of  Eq.  (35)  over  the  new  region 
AS.  We  note  that  J(r,t)  in  this  integral  consists  only  of  tlie  term 
since  the  incident  field  has  just  passed  over  each  patch  in  this  region  so 
that  there  are  no  contributions  from  the  surface  currents  at  other  patches. 

We  can  thus  write  the  first  order  physical  optics  contribution  of  the  new 
region  entirely  in  terms  of  known  quantities  by  rewriting  Eq.  (35). 


where  r H is  the  known  ramp  response  of  the  surface  and  where  r H , and 
o ft-  Q 

r H are  the  contributions  to  the  far  field  due  to  the  known  region  (S-As) . 
o c 

The  left-hand  side  of  Eq.  (41)  can  be  integrated  directly,  as  shown  in  the 
next  subsection.  Since  the  self-term  correction  factor  (as  it  appears  in 
Eq.  37)  is  dependent  upon  the  surface  geometry  and  the  polarization,  we  can 
obtain  the  geometry  characteristics  by  comparison  of  the  responses  at  the 
two  polarizations. 


We  can  see  now  how  the  inversion  process  would  proceed : 

1.  For  the  specular  point  tlie  riglit-hand  side  of  Eq.  (41) 
is  simply  the  leading  edge  of  the  ramp  response,  hence 
solution  of  the  integral  will  yield  the  geometric  prop- 
erties of  the  region  surrounding  the  specular  point. 

2.  {^lowing  tlie  region  (S-As)  , update  the  correction  currents 
J up  to  the  retarded  t^me  corresponding  to  the  new  far 
field  time  t^;  compute  J and  J„„,  at  the  retarded  time 

f r 

and  their  contributions  H and  H to  the  far  field. 

c roi 

3.  .Subtract  these  contri)Mitions  from  tlie  known  ramo  response. 
Do  this  for  both  polarizations. 

4.  From  the  remaining  contributions  to  the  far  field  deduce 
the  geometric  properties  of  the  new  region  As. 
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5.  Tlie  entire  surface  S is  now  known.  Step  to  a later 
time  t^  and  repeat  for  a new  region  /.S . 

It  remains  to  relate  the  characteristics  of  t)ie  surface  to  tlie 
leading  edge  of  its  first  order  ramp  response. 


5.3  FIRST  ORDER  RAMP  Ri:SPONSE  OF  CURVILIUEAR  PATCH 

We  want  to  reformulate  the  following  integral,  which  represents  the 
known  contribution  to  the  far  field  of  the  new  unknown  region  uS: 


I(r  ,t  J 
o f 


_i  _a 

4it  3t 


As 


o 


(42) 


where  AS  is  a region  of  as  yet  unknown  geometry  and  where  t^  is  such  that 

As  is  just  covered  by  the  incident  field  at  the  corresponding  retarded  times 

t = t + a • r.  This  is  illustrated  in  Figure  33.  VJe  need  to  find  an  ex- 
f r 

prossion  for  this  integral  in  terras  of  the  geometric  properties  of  AS. 

By  the  result  of  Appendix  8.3  we  have  that 

> > 1 ^ *-inc  “*■  1 "^inc  -*■ 

J (r,t)  = J (r  ,t)  = ~ (2a  x H (r  ,t)  H - J (r,t)  . 

POl  Y n Y 

- >■ 

To  determine  Y it  in  nece.ssary  to  write  J explicitly  in  its  u,v  components 
(.see  Appendix  8.3); 


J (r,t; 
u 


me  * 

J ( , t ) 
u 


1 - C 


(r,t) 

J (r,t)  = — 1 
V 1 + 1 


(43) 


wnero 


K - K 
^ _ _ii y 


and  wliorc 


(44) 


R is  tlie  radius  of  the  self-natch.  Tlie  critical  point  here  is  that  the 
o 

value  of  R is  proportional  to  the  value  of  the  incident  ramp  at  (r,t)  . Uf' 

^ 1 n o * ^ 

see  this  as  follows:  Tlie  magnitude  of  J at  (r,t)  is  given  by  Kq.  34. 

J^'^'^(r,t)  = (a  • r)  + t,  while  the  magnitude  at  an  integration  point  of  the 
o 
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self-patch  is  the  retarded  excitation: 


inc 


(r",T)  = (a  • r")  + T 
r 

o 


or,  using  the  usual  definitions  R = r-r^  , x = t-R: 


^inc,->-^  . ,inc  ,->■  . ,-s  :X, 

J (r  ,T)  = J (r,t)  - R - (a  • R)  . 

r 

o 


Note  that  if  the  curved  patch  is  perpendicular  to  a^  , the  dot  product  term 
is  small  in  comparison  to  R.  On  the  other  hand,  if  ?he  patch  is  nearly  tan- 
gent  to  a^  , (a^  * R) ~ ± R.  If  we  average  the  effects  of  2 points  on  oppo- 
site sides°of  r,°we  can  thus  simplify  the  above  as 


= j'^’^^(r,t)  - R. 


(45) 


me 


This  is  greater  than  zero  for  R<J  (r,t) , hence  the  radius  of  the  self- 
patch is 


R = J^‘"'^(r,t) 
o 


This  is  illustrated  in  Fig.  34:  The  self-patch  that  contributes  to  the  sur- 

face  current  at  (r,t)  for  the  leading  edge  of  the  incident  ramp  is  only  a 
small  region  about  r.  We  consider  now  two  directions  of  polarization  of 


a and  a , where  a = a x a 

X p p X r 


(a^  and  a^  are  in  the  yz  plane).  In  Appendix  8.4  it  is  sliown  that  if  we 
transform  °J  from  the  u,v  coordinate  system  of  Eqs.  (43)  to  the  x,p  coordi- 
nate system  and  form  the  cross  product  required  in  Eq.  (42),  we  can  write 


.l(r,t)xa  =2(a*a  ) 

r nr 

o o 


me  .N 

/H  (r,t)a 


1 + C 


H^”"'(r,t)a 

X X 


(46) 
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leading  edge  of 

INCIDENT  RAMP 


LEADING  EDGE 
OF  INCIDENT  RAMP  AT 
PREVIOUS  ITERATION  STEP 


BY  PREVIOUS 
ITERATIONS 


FIG.  33  Unknown  region  A S for  body  with  rotational  symmetry. 


EFFECTIVE  SELF 
PATCH  AREA 
ABOUT  POINT? 


FIG.  34  Self  patch  integration  at  point  inside  A S. 
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Hero  certain  crosa-torms  are  neylectod,  which  cancel  when  the  above  is  into- 
qratod  over  a surface  with  plane  symmetry  about  the  yz  plane  (as  is  the  case 
for  our  bodies  of  revolution) . 


Let  us  concentrate  on  the  a component,  and  substitute  Lc}.  (46) 
inc  ^ 

and  the  value  of  U in  the  expression  for  the  far  field  contribution  (using 


t = t^  + a • r) : 
f r 

o 


t^  + 2(a  • r) 

1 ^ f ^ 

I (r  ,t  ) = — 7—  I — 2 (a  • a ) dS 

p o f 271  3t  ( 1 + r ) n r 

As  ' ° 


A ^ r 


I (r  ,t^)  = — ^ — - dS 

' 271  3 1 J (1  + ^)  pro] 


p o f 


(47) 


As 


Here  we  have  written  dS  . = (a  • a ) dS , the  projected  area  and 

ppoj  nr 

t = t + 2(a  • r)  . Note  that  t^  = 0°at  tlie  far  edge  of  AS  where  the  Icad- 

f r 

ing  edge  of  tlie  ramp  has  just  arrived  and  is  >0  in  the  rest  of  .Is.  Also 


note  that 


(K  -K  ) 
u V 

C = ^ • t 


Since  tlie  contribution  to  the  integral  at  the  edge  goes  to  zero,  we  neglect 
the  effect  of  time  variation  of  the  integration  limits  and  taicc  the  time 
derivative  inside  tlie  integral: 


I (r  ,t ) = ^ r 

p o f 2tt  , 2 


(4«) 


As 


(1  + A) 


Tliis  integral  can  bo  performed  over  t , by  writincj  dS  . in  terms  of  t . 

pro] 


The  limits  of  t are  t =0,  corresponding  to  tlie  far  edge  of  AS, 


and  t”  = At  = 2A(a  ' r)  , corresponding  to  tlie  near  edge  of  l.S  (winch  is  the 
o 


i. 
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edge  of  the  known  region).  Expand  dS^^^^  from  a point,  e,  on  the  near  edge: 


dS  . / dS 

prop  ^ 1 __E 
dt  dt 


(At  - t')  + 


(At  - t")  + ... 


and  expand 


= 1 - 2?  + 3? 


(1  + ^)' 


Integrating  and  retaining  only  second  order  terms 


; dS  . d^S  . \ .2  (K  - K ) dS 

2-nr  = i ■■  At  + I ^ y ^ ^ 

^ ^p  \dt  j ^ ,2  / 2 2 dt 

^ \ /e  \dt  'e  ' 


^ At 

7 2 


Similarly,  for  the  other  polarization,  we  obtain 


2TTI  = . 

X \ dt 


.2  (K  -K  ) dS^ 
At  u V £ 

2 2 dt 


The  first  two  terms  of  these  are  the  projected  area  of  the  new  region  As, 


I + I 

X t 


which  is  the  polarization-independent  physical  optics  term.  Also,  approxi- 
mately, subtract  Eqs.  (49)  and  (50)  obtaining 

^ - I /I 

K -K  = . ^ 

V ■ I + I At 

\ X p • ' 


To  obtain  K ,K  we  need  anotlier  relationship.  At  the  specular  point  we  can 
u V 


use  the  fact  that  (see  Appendix  8.2) 


proj  _ TT 

" /K^K  ' 
U V 
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resulting  in  a quadratic  equation  in  K or  when  combined  with  £'q.  (52). 

On  the  preceding  page,  it  is  assumed  that  ^0.  In  Appendix 
8.6  the  case  of  a specular  return  from  a cylindrical  surface  is  described. 
At  a non-specular  point,  we  malte  use  of  our  assumption  that  S is  a surface 
of  revolution  about  the  2-axes.  In  Appendix  8.2  it  is  shown  that 


‘V  = - 


p/l  + 


dp 

dz  . 


(54) 


We  )cnow  now  the  projected  area  of  the  new  region,  from  Eq.  (41),  and  we  )<now 
Its  depth  from  our  choice  of  t^.  Using  the  expression  given  in  Appendix  8.5 
relating  the  projected  area  of  a section  of  a body  of  revolution  to  its  di- 
mensions we  can  solve  for  t)ie  D,z  coordinates  of  the  new  region  for  the  ori- 

2 2 

entation  dO,^dz,  and  hence  for  K and  K . (We  can  also  find  d 0,^dz  from  the 

u V 

relation 


K 

u 


1 + 


do 

dz 


3/; 


(55) 


which  is  useful  for  tlie  initial  estimate  of  the  location  of  the  next  un)tnown 
region . ) 

Tlie  new  region  is  thus  completely  specified  by  this  analysis.  Some 
com-ments  are  required  here  indicating  the  assumptions  made  in  this  process; 
the  new  region  consists  of  two  parts,  one  to  each  side  of  the  l^ncwn  region. 
Only  for  broadside  or  axial  incidence  will  these  two  regions  have  the  same 
geometric  properties.  For  general  oblique  incidence,  or  for  a body  wit;iout 
the  left-right  s^rnimetry  assumed  here,  it  would  be  required  to  use  a somew'nat 
more  complicated  process  of  curve  fitting  to  find  the  best  extension  of  fne 
contour  in  each  direction  consistent  wit;i  the  calculated  average  values  of 
the  surface  properties.  Such  a process  appears  to  be  a natural  extension  of 
the  methods  presented  iwre,  althougli  it  'nas  not  been  implemented. 
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5.4  NU.'IKRICAL  SOLUTION  P.ROCLDUI-IE 


i 


i 

t 


[ 


By  way  of  summary,  we  outline  here  the  stop  by  step  inversion  pro- 
cedure, as  implemented  in  the  computer  program  INVRS. 


1.  Input  of  the  ramp  responses  H (t)  and  H (t) . These  are 
obtained  from  the  space-time  integral^^  equation  solu- 
tions, as  described  elsewhere  (4),  augmented  by  the  ideal 
specular  response  at  t = 0 (see  Appendix  8.6),  and  given 
numerically  at  regular  intervals  of  t. 

2 . Analysis  q£  Specular  Return 

Equations  (41)  througli  (45)  are  applied  at  tlie  leading 
edge  of  the  response  in  the  limit  as  /:t  *0  (i.e.  we 
actually  calculate  the  derivatives  of  the  responses) , to 
establish  the  curvatures  at  the  specular  point.  The 
accuracy  of  this  process  is  limited  by  the  fineness  with 
whicli  the  ramp  response  is  specified.  On  the  assumption 
that  the  target  is  a body  of  revolution,  and  Icnowing  the 
angle  of  incidence,  we  can  establish  bands  along  the  ra- 
dius of  revolution  tlirough  this  point.  This  is  our  first 
Jenown  region. 

3.  Extrapolate  from  the  )criOwn  region  by  a small  distance 

to  each  side  to  obtain  the  approximate  center  of  the  next 
bands.  The  part  of  these  bands  that  will  be  illuminated 
in  accordance  with  step  (4)  is  our  new  region  Ls. 


4.  Choose  the  far  field  time,  t,,,  such  tliat  the  leading  edge 
of  the  incident  ramp  just  reaches  the  far  edge  in  the  y-z 
plane  of  the  new  region  AS  at  the  retarded  time 


t 


t + a 
f r 


'edge ' 


5.  Update  the  surface  correction  currents  at  all  points,  i,  in 
the  Icnovm  region  up  to  t'.ieir  respective  retarded  times 

t.  = t..  + a • r,  , using  Eq.  (39). 

1 f r 1 ' 

o 

6.  Calculate  the  contributions  of  the  V;nown  region  to  t’ne  far 
field  at  time  t^  (for  cacli  polarization)  using  Eqs.  (39) 
and  (40)  . Subtract  these  from  the  !<nown  rcsp<3nses;  t)>e 
remainder  is  due  to  tlie  new  region. 

7.  Apply  the  expressions  (41),  (42),  (44)  and  tlic  result  of 

Aptjcndix  8.5,  to  characterize  t!ie  geometry  of  t'no  .low  region. 


J 
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now  extended  our  ki'iown  recjion  and  continue  Erom  point 
(3)  until  the  entire  surface  is  covered.  (For  a body 
or  revolution  it  is  only  necessary  to  pass  halfway  over 
the  target.) 

At  point  (2),  the  analysis  of  the  specular  return,  the  determina- 
tion is  made  of  the  nature  of  tiie  specular  surface.  That  is,  whether  one  or 
both  of  the  curvatures  are  zero  (See  Appendix  8.6).  In  the  case  of  the  sphere- 

capped  cylinder  viewed  on  broadside,  for  exa.mplc  one  of  the  curvatures  is 

zero,  and  tlie  initial  known  region  is  the  entire  cylinder  body. 

The  core  of  the  process,  point  (5) , the  updating  of  tiie  surface 
currents  is  similar  to  the  usual  space-time  integral  equation  scattering 

solution  [4,5];  except  that  here  tlie  time-derivative  of  the  correction  cur- 

rents are  calculated  rather  tlian  the  usual  total  surface  current.  Also,  the 
correction  currents  are  calculated  at  the  jxact  retarded  time  corresponding 
to  the  far  field  time  under  consideration,  in  addition  to  the  calculations 
at  the  regular  time  grid.  The  patch  interaction  of  the  correction  currents 
(Eq.  39)  are  calculated  at  tiie  retarded  times  T by  interpolation  on  the  time 
grid,  whereas  t.he  effects  of  the  incident  field  are  calculated  exactly. 


5.5  INVERSE  SCATTERING  RESULTS 

We  present  here  the  results  of  exercising  the  computer  program 
INVRS  on  the  test  cases  of  a sphere,  a prolate  spheroid  and  a sphere-capped 
cylinder,  all  viewed  at  broadside. 

A quadrant  of  the  cross-section  of  a sphere  is  compared  with  tlio 
result  from  the  inverse  scattering  calculation  in  Fig.  35.  The  arrow  indi- 
cates tiio  direction  of  incidence.  Crosses  indicate  tlie  sample  points  of  tlie 
calculated  contour.  The  agreement  is  very  good  up  to  tlie  last  jioint  calcu- 
lated near  the  axis  of  symmetry.  At  this  point,  the  projected  area  of  tlie 
new  region  is  much  smaller  than  its  surface  area;  hence  a small  error  in  the 
calculation  of  the  correction  currents  loads  to  a larger  erroi  in  the  surface 
contour.  A perspective  drav>;ing  is  given  in  Fig.  36. 
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FIG.  35  Inverse  scattering  results  — cross-section  of  sphere  quadrant. 


DIRECTION  OF 
INCIDENCE 


(a)  True  Contour 


(h)  Calculated  Contour 


FIG.  36  Inverse  scattering  results  — sphere. 
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A similar  set  of  drawings  is  pre---ente(J  in  I'lg.,.  37  and  38  for  a 
prolate  splieroid  of  2:1  axis  ratio  with  its  major  axis  along  the  axis  of 
rotation.  T!ie  contour  acjroemciit  for  tliis  case  is  very  good,  altiiougli  there 
is  a 5%  overestimate  of  ttie  size.  Tlii.s  small  error  is  due  to  inaccuracies 
in  the  specular  point  analysi:;,  caused  by  coarseness  in  tlic  numerical  repre- 
sentation of  tile  ramp  ri'sponse. 

Tiio  results  for  the  sphere-capped  cylinder  are  represented  in 
’’iqs.  39  and  40.  The  cylinder  length  was  3 times  its  diameter.  The  results 
of  this  case  were  less  satisfactory  than  for  the  other  test  cases.  There 
are  several  prolilems  in  the  analysis  of  tills  case:  The  analysis  of  the 

specular  point  (as  descr  ibc-ii  in  .\ppendi.x  8.6)  requires  sep/arating  the  nart 
of  the  response  due  to  the  cylinder  (which  has  a dependence  of  half-integral 
powers  of  time)  from  the  response  due  to  the  smoothly  curved  cap  (which  has 
a dependence  of  integral  power;!  of  time) . The  numerical  representation  of 
the  ramp  response  available  was  too  coarse  to  permit  an  accurate  analysis 
of  this  type.  Instead,  reasonable  values  for  the  dimensions  of  the  cylin- 
der as  might  be  expected  from  the  analysis  of  more  accurate  response  data, 
were  utilized  in  the  solution.  The  remainder  of  the  contour  generation  is 
reasonably  good  until  the  region  closest  to  the  axis  of  symmetry  is  reached. 
Here  we  have  the  same  effect  noted  for  the  case  of  the  sohere,  namely  tliat 
small  inaccuracies  in  the  calculation  of  t)ie  correction  currents  can  pro- 
duce large  errors  in  that  part  of  the  contour  v;hich  is  nearly  parallel  to 
the  direction  of  observation.  The  error  is  larger  in  the  case  of  t.he  cyl- 
inder tlian  for  the  sphere,  because  of  the  much  larger  surface  area  over 
v/hich  the  surface  current  error  is  accumulated. 

In  summary:  It  has  been  demonstrated  tliat  inverse  scattering  of 

rotationally  symmetric  bodies  using  the  space-time  integral  equation  an- 
proach  yields  reasonably  accurate  results  for  oblique  incidence.  Accurate 
specular  analysis  requires  the  representation  of  the  response  at  intervals 
of  At = .05*R,  where  R is  the  body  radius  (At=0.2"R  was  available  tor  tnese 
test  cases) . Inaccuracies  in  the  surface  correction  currents  can  j'robably 
be  improved  by  greater  attention  to  the  details  of  surface?  excitation  at 
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FIG.  37  Inverse  scattering  results  — cross-section  of  prolate  spheroid. 


(b)  Calculated  Contour 


FIG,  38  Inverse  scattering  results  - prolate  spheroid.  | 
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FIG.  39  Inverse  scattering  results  - cross-section  of  sphere-capped  cylinder. 
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tile  instant  of  the  t amj  initial  inciiience  as  was  lioni-  loi  thi.-  ea;e  of  axial 
ineiilonce  in  Section  4.  Also,  although  tiie  test,  cases  were  for  broad.si'le 
i nc ivieticf' , it  has  bei'ii  inilicatod  in  the  Ixjdy  of  the  report  that  tlie  method 
can  be  extended  to  arbitrary,  but  knov;n,  anylo  of  incidence  on  surfaces  of 


rotation. 


■H 
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SUHI'JVHY  AND  CON’CLUSIOtlS 

In  this  study  the  technique  for  comiJUtinq  the  smoothed  impulse  re- 
sponse iias  been  extended  to  t)ie  case  of  scattering  from  simjJle  aircraft 
models.  Those  models  consist  of  a cylindrical  fuselage  and  rectangular 
plate  wings,  stabilizers,  and  rudder.  To  the  knowledge  of  these  writers 
this  is  the  first  time  that  a successful  solution  for  this  problem  has  been 
reported.  Scconilly,  direct  time  domain  techniques  have  been  developed  for 
the  solution  of  the  inverse  scattering  problem  involving  rotationally  sym- 
metric targets.  Bath  the  case  of  a.xial  incidence  and  the  case  of  oblique 
incidence  were  considered.  The  use  of  polarization  information  enabled  tlie 
develop'ment  of  a solution  for  the  case  of  oblique  incidence.  Finally,  tech- 
niques for  computing  the  impulse  response  have  been  developed  for  the  case 
of  cylinders  with  fins  attached  with  non-axial  incidence.  The  foundation 
of  these  techniques  is  tiio  si>ace-time  integral  < juatio:i  .ipnroach  to  the 
solution  of  the  electromagnetic  scattering  problem.  This  coupled  with  the 
impulse  response  augmentation  technique’,  is  used  to  yield  tlie  total  impulse 
response  or  the  frequency  response  over  the  entire  spectrum.  This,  in  turn, 
can  be  used  to  compute  the  response  of  a target  due  to  any  incide.nt  radar 
waveform,  regardless  of  wave  shape  or  carrier  frequency. 

Numerical  procedures  v^rere  developed  for  tlie  solution  of  a space- 
time  integral  equation  for  computing  the  smoothed  impulse  response  of  simple 
aircraft  models  consisting  of  a cylindrical  fuselage  with  rectangular  plate 
wings,  stabilizers,  and  rudder.  These  procedures  were  extensions  of  pre- 
viously developed  techniques  for  computing  the  smoothed  impulse  response 
of  a cylinder  with  single  fins  attached.  They  consist  of  the  sinultanoou.s 
solution  of  two  space-time  integral  equations  for  the  surface  current.  One, 
an  H-fiold  type  equation,  is  for  the  cylindrical  fuselage  surface  current 
and  it  includc.s  t!io  effect  of  currents  on  the  fuselage  and  currents  on  the 
wings,  stabilizer,  and  rudder.  The  second  equation,  an  F-ficld  type  equa- 
tion, is  for  the  current  on  tlie  wings,  stabilizer,  and  rudder  and  it  in- 
cludes the  effect  of  both  the  wing,  stabilizer  and  rudder  cut  rents  and  tlie 
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iu;;<'Unio  cui't  t-ntii . Finally,  at  tliu  bcjundary  between  the  tins  an  i the  fune- 
laeje  the  phy.sical  Ijoundary  conditions  are  ajijjlied  to  specify  tiie  current  in 
that  teqion.  Thin  technique  was  verified  l>y  direct  tiiiK,  lioinain  calculations 
to2’  a nui:ib(.'r  of  simple'  canes  and  found  to  be  in  qood  aqre't.-nent . The  tech- 
nique wa.s  demonstrati'd  on  an  aircraft  model  that  consisted  of  a cyliiuirical 
iuselaqe  with  rectanqular  fins  representiny  winys,  stabilizers  and  rudder. 
These  results  are  particularly  encourayiny  and  the  tecliniqueu  developed  here 
olioulci  bo  t'xtondetl  to  more  refined  aircraft  models.  T!u?se  more  refined 
racviels  could  consist  of  smooth  fuselayc  representations  which  more  closely 
approximate  the  real  aircraft.  In  addition  tlie  wings,  stabilizers,  and  rud- 
ders could  be  represented  by  quadrilateral  flat  plates  instead  of  the  now 
liraitiny  rectanqular  flat  plate  yeometrios. 

The  impulse  response  auymentation  teclini(}ue  was  successfully  ex- 
tended to  the  case  of  scatteriny  from  a cylinder  witii  a single  £>air  of  fins 
attached  for  oblique  incidence.  The  augmentation  functions  for  this  target 
were  based  on  the  results  previously  obtained  for  the  cases  of  square  plates 
and  riylit  circular  cylinder.s.  By  using  a combination  of  these  two  sets  of 
au>;r.nentation  functions  the  necessary  augmentation  fund  ion  for  the  case  of 
a cylinder  with  fins  attached  was  obtained.  Results  v;ere  computed  for  a 
cylinder  with  a length  to  diameter  ratio  of  2:1  and  with  one  set  of  square 
fins  whose  size  was  equal  to  the  cylinder  diameter  attaclied.  These  results 
wore  obtained  for  aspect  angles  of  0^,  30°,  60°,  90°,  120°,  150°,  and  180°. 
The  resulting  impulse  responses  clearly  show  tlic  singular  contributions  from 
the  various  scattering  centers  on  this  more  complex  target.  The  cc'mpanion 
frequency  responses  display  ripples  that  correspond  to  the  interference 
between  tiie  returns  from  these  different  scattering  centers.  The  results 
of  this  technique  continue  to  be  encouraging  and  -hould  bo  extended  to  more 
co.'f’.plex  targets  and  to  non-rcctangular  fin  plate  targets. 

A time  (.iomain  approach  to  the  inverse  scattering  problem  was  de- 
veloped and  demonstrated  in  this  study.  This  apjiroacli  starts  with  a space 
time  integral  equation  that  represents  the  inversion  process.  Previously, 
an  iterative  techniejue  was  developed  for  tlie  solution  of  this  inversion 
equation  and  applied  to  four  tarejet  geometry  c.ises.  In  this  effort  a direct 
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t iwo  ilomain  solution  ti’chnitjuo  was  developed  for  the  solution  of  this  inver- 
sion ejuation  ami  applied  to  the  same  four  targets  that  were  examined  pre- 
vious!'/. Results  wore  obtained  for  tlie  case  of  axial  incidence  for  a spiiere, 
spiiere-cap  c/linder,  and  a flat-end  sphere-cap  cylinder.  The  direct  solution 
technique  results  compared  wi'll  with  tiiose  obtained  using  the  iterative  solu- 
tion technique  and  were  in  close  agreement  to  the  actual  target  geometries. 
However,  tlie  major  advantage  of  the  direct  solution  tecliniquo  is  that  it 
requires  substantially  less  computer  time  to  obtain  the  .solution.  In  addi- 
tion to  the  case  of  axial  incidence,  the  case  of  oblique  incidence  was  also 
considered  for  the  inverse  problem.  In  the  case  of  oblique  incidence  the 
direct  time  domain  techniques  were  developed  for  the  solution  of  the  inverse 
scattering  problem  using  polarization  information  together  with  the  target's 
ramp  response.  As  an  example  of  tliese  tech.niques,  the  case  of  broadside 
incidence  on  three  targets  was  considered.  The  results  were  obtained  for 
the  case  of  a sphere,  a sphere-cap  cylinder,  and  a prolate  spheroid.  These 
results  yielded  a geometry  estimate  which  was  in  close  agreement  witli  the 
actual  geometry  except  at  the  ends  of  the  targets.  It  was  found  that  this 
technique  requires  only  a single  observation  angle  for  the  case  of  broadside 
incidence . 
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SKCTIOn  8 
APFC’.IDIX 


B.l  IMPROVED  SPACE-TlMi;  INTEGRAL  EQUATION  SOLUTION  TECHNIQUE 


To  implement  the  solution  for  the  inverse  scattering  problem  des- 
cribed in  Section  4,  it  is  necessary  to  compute  tlie  correction  currents, 

J , for  the  surface.  T.heso  currents  represent  the  effect  at  a given  point 
CR 

on  the  surface  of  currents  flov/ing  at  other  surface  points.  In  this  compu- 
tation use  is  made  of  the  rotational  symmetry  of  the  target.  This  allows 
for  a substantial  reduction  in  tiie  computer  time  and  memory  requirements. 


The  geometry  of  the  rotational ly  symmetric  scattering  problem  is 
displayed  in  Fig.  23.  For  convenience  the  polarization  of  tlie  incident 
field  was  taken  to  be  vertical,  since  the  (p  reference  is  arbitrary  in  rota- 
tionally  symmetric  problems.  The  scatterer  surface  is  defined  by  the  con- 
tour 

P = p(z) 


where  p (the  usual  cylindrical  coordinate)  is  the  distance  from  the  z-axis. 
The  contours,  produced  in  the  y = 0 plane,  starts  at  the  most  positive  point 
in  2 and  moves  along  the  curve  in  the  negative  z-direction.  The  surface  is 
then  generated  by  rotating  the  contour  about  the  r-.-axis.  The  unit  vectors 
of  interest  are 

a = the  unit  vector  normal  to  the  surface 
n 

a = the  unit  vector  tangent  to  the  surface  which  lies  in  the 
olane  generated  by  t)u'  z-axis  and  tlie  position  vector 
r.  (This  corresponds  to  a-^  in  the  spherical  coordinate 
system. ) 

a = the  unit  v’ector  tangent  to  the  surface  which  is  perpen- 
^ dicular  to  a^.  (This  corresponds  to  a in  the  spherical 
coordinate  system.) 

It  is  also  convenient  to  define 
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g 


1 


1 + 


/dp\' 
V dz/' 


dz 


^ ' 1^) 

dz  .1 


which  yield  the  following  relations 


= - g a - ga 
2 p z 


3 4> 


a = ga  - g a 
n p p z 


The  expression  for  the  total  surface  current  J , which  was  derived 
in  Section  4.1,  is  repeated  here  for  the  case  of  an  incident  ramp  waveform 


J„(r,t)  = 2a  X H^(r,t)  + j (r,t) 
R n R CR 


(56) 


where 


1 

27T 


a X 
n 


i*iif  ’‘*R 


R 


T = t - R 


By  virtue  of  the  rotational  symmetry  in  this  problem,  the  total  surface  cur- 
rent can  also  be  written  as 


J (r,t)  = a.,  J 01  ,t)  cos  4>  + a,  (2.  ,t)  sin  (J)  (57) 

R 2 R^  z J z 


where  2 is  the  position  coordinate  in  the  a direction.  Combining  these 

two  equations  and  solving  for  J yields  the  coupled  integral  equations: 

CR 
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2-rr 


J ^a 

' J.  + -1  -l" 

3 2 3t 

S 

R R 

ds' 

^2  ^ 


T = t - R 


r f 

"-  + -^1 

J b 

„3  2 3t 

S 

lR  R .1 

T = t - R 


(58a) 


_J^ 

2tt 


3 2 9t  ' “^R 

""  : r"^  R^  I ^2  ^ 


T = t - R 


r 


1 3 i 

! R~  R^  ' ' ; "3 


+ ~ f f , ! ^ -5^  i J (il",!)  dS' 

2tt  J d ' 3 _2  3t  . R,  z 


(58b) 


where 


T = t - R 


f = - Ijg^p^  + g'  (z-z^))  cos  4>'  - pg^  i cos 
a p ; j 


^ 2 ^ 
= - (z  - z ) sin  (j) 


f = l-ggp  + ggp-gg  (z-z  ) sin 


2.  - 


f = I gp  - g (z-z^)lcos  <j)  - gp'  cos  4>'  • 

alp  J 

The  initial  approach  ta)cen  for  the  solution  of  Eq.  (58)  was  that 
employed  in  the  direct  time  domain  solution.  A grid  is  structured,  as  des- 
cribed in  Section  4.2,  on  the  portxon  of  the  target  whoso  geometry  has  been 
determined.  The  incident  field  is  taken  to  bo  a ramp  defined  by 


H^(z,t)  = { 


0 ; t + z ^ 0 

I t+z  ; t+z>0 
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For  the  numerical  calculations  Eq.  (58)  is  represented  as 


J „ (z. ,t) 
CR^  1 


2tt 


j ^ 


1 X 9 

9xf  R^  1 ^ 


f . 


ij£  ijX. 


. J.  V <r  f 1 1 

2tt  / L bij£  3 ■"  2 

j SL  t^jJl 


9 I 


(t.,T)  AS^ 

ijil  j 


J (z . , t) 
CR^  1 


= + 


5^)'V 


j 


f • • 


^ 'J.,  (^.,r)  hs. 


cij*-  ^,o3  " „2  9xf"R^'"'j 


R.  . . R.  . 0 
^ 13J6  iji 


T = t - R 


ija 


7 


j i 


; ijil  ijH  J 

r = t - R . 


ij£ 


where 


z.  designates  the  band  on  whicli  the  observer  is  located. 


z.  designates  the  band  on  wliicli  the  integration  point  is  located. 


designates  the  patch  on  band  j on  which  the  integration  point 
is  located. 


^ai  j X- 


4>o  < N p-  + g . (z.  - z.)  I cos  !}>.  - p.g.  > 
I L 3 3 P3  1 3 J ^ 13/ 


^bi  j 1 


sin  d>„  (z  . - z . ) 
I 1 J 


2 

f . . , = sin  4>olg  ,g.p.  -g  .g.p.  -g  .g  . (z.-z.) 
ciji  I P3  1 1 pi  3 3 pi  P3  1 3 - 


^dijii  = 


R|.£  = the  distance  between  the  observation  fKDint  and  the  inte- 

3 gration  jxjint. 


(59a) 


(59b) 
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= the  area  of  patch  I . 


The  time  differentiation  and  interpolation  necessary  for  the  eval- 
uation of  the  integrands  in  Eq.  (59)  are  performed  numerically  by  represent- 
ing the  surface  current  with  a fourtli-order  polynomial.  In  order  to  achieve 
the  best  accuracy,  the  five  points  used  for  the  representation  were  chosen 
such  that  the  current  would  be  evaluated  as  near  as  possible  to  the  middle 

of  them.  The  correction  currents  J (z.,t)  are  computed  for  the  time  at 

CR  i 

which  the  leading  edge  of  the  incident  field  reaches  z^.  This  insures  that 

all  non-zero  contributions  to  J come  from  portions  of  the  target  whose 

CR 

geometry  is  already  known. 

For  the  calculation  of  p(z,t)  it  is  necessary  to  evaluate  the  time 
derivative  of  integral 

J N 

I(r,t)  = j*  {J  (r',T)  X a '•  ds'  (60) 

,1  i CR  IT 


For  the  numerical  implementation  I(r,t)  is  represented  as 
I(z,t)  = y ^ la.j^  ^ 


where 


^i£  “ • 


The  time  derivative  of  the  sum  l(z,t)  is  calculated  using  the  central  deri- 


vative formula 


3l (z,t)  ^ I (z,t  + c)  - I (z,t  - e) 
at  2e 


whore  e = .01. 
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The  results  obtained  using  this  numerical  implementation  gave 
solutions  which  were  quite  close  to  the  actual  target  contour.  They  did 
not,  however,  reproduce  the  contour  exactly.  In  an  attempt  to  improve  the 
results,  it  was  found  that  a more  precise  computational  method  was  needed 
to  account  for  the  surface  current  interaction  caused  by  the  leading  edge 
of  the  incident  ramp  waveform.  The  effect  of  the  derivative  in  the  expres- 
sion for  J combined  with  the  derivative  necessary  for  the  computation 
CR 

of  p(z,t)  is  to  necessitate  integrating  an  impulse  in  discrete  space-time 
steps.  This  discontinuity  occurs  only  at  the  leading  edge  of  the  incident 
ramp  since  the  second  derivative  of  the  ramp  at  (z,t)  is  zero  after  the 
leading  edge  passes. 


derivative 


To  obtain  a new  numerical  representation  the  expression  for  the 

of  the  correction  currents,  3j  /Dt  was  expanded,  to  yield 

C K 


3t 


Jl 

2tt  3t 


i 


- — ' j 
R 3t  i PO 


(r  ,T) 


+ J (r 
CR 


X 


a 


T = t - R 


(62) 


where 


2a  X (r^ ,T) . 
n R 


Since  the  discontinuity  occurs  only  at  the  leading  edge,  the  contribution 
to  '^PQ  only  portion  of  the  equation  which  must  bo 

refined.  It  is  necessary,  then,  to  obtain  a numerical  representation  for 


3t 


J.J.  r 

2TT  3t  j 

s 


a 

n 


1 _3 
R 3t  , 


(r ', 1 ) 


X = t - R 


(63) 


For  the  numerical  model  the  target  is  represented  by  a grid  of 
circular  patches.  Each  patch  is  considered  to  bo  in  one  of  ttiree  states  at 
time  T depending  on  whether  the  leading  edge  of  the  incident  field  tins  not 
yet  reached,  is  located  in,  or  has  passed  tlirough  tlic  patch.  If  tlie 
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incident  field  has  not  yet  reached  the  patch,  then  J and,  consequently, 

-*■ 

^V3t  are  zero.  If  the  incident  field  has  passed  through  the  patch,  then 

2 

there  is  a contribution  only  from  the  1/R  term  in  Eq.  (63)  since  the 
3 Jp^/Bt  is  zero.  For  this  case,  Eq.  (63)  can  be  represented  numerically 
in  terms  of  its  two  components  as 


j z 


(f  . . „ + f,  . . „ 
aijZ  hijZ 


^p.’  1" 

^ '^ijA 


3j  (z.,t) 

CRl^  1 


^^cijJl  ^dij£ 


“^p.*  3 ^^Z 


If  the  leading  edge  of  the  incident  field  is  within  tlie  patch, 
then  Eq.  (63)  can  be  expressed  in  terms  of  the  fraction  of  the  total  patch 
area  which  is  encompassed  by  the  incident  field.  For  this  case,  the  com- 
ponents of  3j  /3t  can  be  written  as 
CRi. 


1 V V ^ , r 1 3a  1 3a"^  , 

Ti  /.  ^ ^aijJL  ^bij«,  ' '^p.’  3 3t  2 3t  ■ 


1 V V ,r  V f 1 3a  1 3a" 

7T  ) ) ^^cijil  ^dijX,  • “^p.^  3 3t  ^ 2 3t  j 

j £ ^ I ij£  ij£  ' 


i /*2  2'  2 -1/  X \ TTa^  dx 

= . X / a - X + a sin  (-j  - — — 


3a"  ^ '2  2 dx 

ar  s 
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a = — 

TT 


I 

I 
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'.R(a  *a  ) - t(a  *a  ) - /[R(a  *a  ) - t(a  *a  - fl  - (a  -a  )^Tr^  - 

LPR  PPj\^Lpr  PPj  L PPlL  j 


1 - (a^^  • a )• 
P P 


dx 

dt 


(a^*a_)  + fR(a  *a^)(a  •a)-t'/,R(a  *a^)-t(a  ’a)"!^-  ’^l-(a  *a_)^  R^-t^' 

. vL  " 


p p 


P R p p 


p R P P 


P P 


1 - (a^  • a ) 

P P 

a = unit  vector  from  the  integration  point  to  the  observation  point 
R 


a^  = unit  vector  in  the  direction  of  propagation  of  the  incident  field 


(a  -a)-a(a  *a  - a *a) 
£ _ R p n p R n p 

P •^v  /V./N  /S 

(a^  - a)-a(a  *a„-a  *a) 

R p n n R n p 


In  the  refinement  of  the  numerical  representation  then,  two  dif- 

ferent  techniques  are  employed.  The  contribution  to  Ppg 

computed  using  Eqs.  (64)  or  (65)  depending  on  the  state  of  the  integration 

- ► -► 

patch.  For  the  contribution  of  J to  /9t  in  Eq.  (62) , the  techniques 

CR  CR 

devised  for  the  direct  approach  are  still  used. 


This  second  approach  succeeded  in  further  improving  the  approxi- 
mations to  the  target  contour.  Although  the  solutions  still  did  not  exactly 
reproduce  the  target  geometry,  they  have  two  basic  advantages  over  the  pre- 
vious approach.  First  the  contours  produced  tend  to  be  more  smooth  than 
those  previously  obtained  and  secondly,  the  techniques  itself  yields  a more 
stable  solution.  The  results  presented  in  Section  4.3  were  obtained  using 
this  second  approach. 


8.2  USEFUL  DIFFERENTIAL  GEOMETRY  REIATIONSHIPS 

Consider  a surface  of  revolution  about  the  z-axis,  where  we  give 
the  radius  of  rotation  p = p(z).  Denote  derivatos  by  p'  = dp/dz,  etc.  Note 
that  p'  lias  the  opposite  sign  of  z and  that  p**  is  always  negative.  A ixjint 
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on  ; ho  sui  t\ioo  IS 


r = (p  cos  if'  , p sin  i+i  , z) 

Cliooso  tho  pr  incipal  coorrlinatcs  u = -z  and  v = 4’-  We  ttien  have 

(-p  cos  4’  » -P'  sin  4’  < "1)  » where  K | r^J  =]+».' 

(-P  sin  4>  . P cos  4>  . 0)  , where  G H 1 1 ^ 

a = r //k  ; a = r / /g  ; 

u u V V 

and  where  an  clement  of  Icnyth  is  qiven  by 

{di.)^  = i:  du^  + G dv^  . 


^r 

hi  ^'u 
iv  ^ V 

and 


Tile  unit  normal  is 


^ ^ cos  iji  , sin  , -p 

a = a X a 

n u V / - 


The  second  order  relations  are 
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w 


Tho  principal  curvatures  are  then 


K = 
u 


K = 

V 


botli  curvatures  are  0 . For  a cylinder  K =0  and 

u 

principal  coordinates  assures  tliat 


i’  r • r = 0 
u V 


and 


M r "a  = 0. 
uv  n 


For  a prolate  spheroii.1 

K ■-  0 . Our  choicH'  of 
v 


- = -l/(p,'F  ) 


i'.i  IHTFOKATION  OVKR  SFLF-PATCH 

The  follciwinq  differs  from  previous  formulations  [4)  by  the  inclu- 
sion of  the  effect  of  retardation  of  tho  surface  cui'rent  over  tho  self-patch. 
The  ('xi  ression  for  tl\o  surface  current  is 

J(l^t)  = -^  I a X 'i-  J(r',T)  X a } ds'  . (Gb) 

271  .1  n R 

S 

Wo  wish  to  evaluate  the  part  of  this  inti'qral  due  to  patch  S^.  containinq  r 
b\'  direct  integration. 

J (r,t)  = IT'  I ax  if  d(r*,7)  x a'  i dS  , (G7) 

C . 2tt  , n R 

S 

e 


vdiero 


.a  i .s  the  unit  normal  at  r 

ti 
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T = t - R 


„ 118 

^ ~ 2 R 3t 

R 


Write  R as  a Taylor  expansion  about  r: 


^ ->•  1 2 1^  2 ^ 

-R=r-r  = r du  + r dv+—  r du  +—  r dv  + r dudv  , (68) 

u V 2 uu  2 vv  uv 

where  u,v  are  principal  coordinates.  Using  the  relationships  of  Appendix 
8.2  we  form; 


.s  1 2 2 

a * a = - — — (EK  du  + GK  dv  ) 
n R 2R  u V 


(69) 


—r  — » ^ 

Separate  J(r  ,i)  into  its  two  components 


J(r,r)=Ja  +Ja  , where 

U U V V 


r"  = r + r du+r  dv,  hence 
u u uu  uv 


r r 

a = a + du  + dv  , and  similarly 

^ ^ /e  /e 


(70) 


r r 

-V  w . uv  . 

a = a + dv  + du  . 

/G 


/g 


Expanding  the  vector  triple  product  in  (67)  and  considering  only  the 
component  to  second  order : 
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a xij‘j'a'xa^'=£j'H(a  •a)a^'[-|(a’a")a  ] 
n y u u R '*■  u (^L  n R u J L n u R j 


C a , EK  du  A 

= ^ (EK  du  + GK  dv  ) ' + a — ~ — + . . . ’ 

u 2R  u V j ^ u R j ■ 


(71) 


2 2 

£j'  (EK  du  - GK  dv  ) . 

u 2R  u V 


In  the  above,  terms  in  odd  powers  of  du  or  dv  have  been  omitted,  because 
their  contributions  will  cancel  when  integrating.  Similarly,  for  the 
component , obtain 


^22 

a X /jj'a'  X a > = (GK  dv  - EK  du  ) . 

n ; vv  R,  v2R  v u 


(72) 


We  can  expand 


3j  (r,t) 


j (r',T)  = J (r,t)  + R 
u u 


9r 


+ (X  - t) 


3j  (r,t) 
u 

3t 


We  now  assume  that  J is  the  same  everywhere  on  the  patch;  ie. 

u 


3j  3j  (r",T)  3j  (r,t) 

u „ 1 u u 

_ = 0 and  — = — 3^ 


then 


J (r',T)  3d  (r',x)  J (r,t) 


r.j  = 

u 


R3t 


Similarly, 


ji  d = 

V 


d (r,t) 

V 


The  quantities  j',  j'  and  a , a can  thus  be  taken  outside  the  integral  of 
u V u V 

Eq.  (67) , which  can  be  rewritten 


1 
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f n 


(r,t)  = — , J (r,t)  a - J (r,t)  a 
47T  1 U U V V 

L 


I 


{ EK 


2 2^ 
du  - GK  dv 

. Y L 


dS' 


R 


(73) 


We  will  integrate  over  a circular  patch  of  radius  R^,  by  letting 
/ e'  du  = R cos  0 ; ^0  dv  = R sin  0 ; and  dS  = RdRd0 


This  integration  yields 


r( 


J (r,t)  = a J (r,t)  - a J (r,t); 

C - u U V V 


K - K (1 

I R 

4 Jo 


(74) 


We  will  have  occasion  to  use  the  abbreviation 


1 K - K . 

' U V 

X,  = — ; R . 

4 o 


(75) 


Rewriting  Eq.  (66) , we  can  substitute  the  above  result: 


J(r,t)  = J£.(r,t)  + — j . . . dS 

S - S_ 


= acj-scj+^  r 

u u V V 2tt  1 


. dS 


S - S, 


where  the  integrals  are  now  over  the  reduced  area,  excluding  the  self-patcli. 
Hence : 


^ 1 1 r 

J 


. . dS 


S - S, 


. 1 1 r 

I 


S - s. 


(76) 


...  dS' 
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It  will  soinot  imes  bt‘  convenient  to  v/rite  this  as 


J(r,t)  = 


...  as  , 


where  , is  undetstood  to  be  evaluated  in  the  above  manner  by  separati.ng  into 
cuiiU'onent:; . (If  one  of  the  component.s  is  zertj,  y ==  1 * i-) 

8.4  TitAUSFOHflATlON  OF  SELF-l’ATO'd  CORRECTION  TO  THE  INCIDENT  H COORDINATE 
■SYSTEM 

The  quantity  to  be  evaluated  in  the  far  field  integral  Kq.  (30) 


J ( r , t ) X a 


We  arc  interested  in  simplifying  the  first  order  physical  optics  contribu- 


'■inc 

J - . -fine  diric 

X a , where  J = 2a  x H 

Y r n 

o 


Y is  defined  in  Appendix  8.3,  and 


a is  the  direction  of  the  far  field  which  is 

If 

o assumed  to  lie  in  the  yz  plane  at  angle  a 
v/ith  the  z axis. 


The  incident  field  H travels  in  the  negative  a direction. 

r 

o 

;^inc  , ,,inc  , , . . ^ 

H (r , t)  = H (r ,t)  • a 


whore 


a = a X a 
p X r 
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I 


and  whoro  either 


a = a ; TM  mode,  or 

H X 


a = a ; TK  mode. 

U p 


,inc.^  ,inc.v 

►me  J a J a 

J u u ^ V V 

Y “ 1 - C 1 + ’. 


Cotisider  the  TM  mode:  a = a 

H X 


me  ^ , ,.ine 

J = 2a  • (a  X a ) H 
u u n X 


^me  ^ „tne 

J = 2a  • (a  X a ) H 

V V n X 


It  is  eonvenient  to  write  out  the  various  unit  veetors  in  Cartesian  form: 


a = (1,  0,  0) 

X 


a = (0,  sm  a,  cos  a) 
r 

o 


a = (0,  - eos  a,  sin  CX) 
P 


a = (n  , n , n ) 
n X y e 


a = (u  , u , u ) 
u X y r; 


a = (v  , V , 0)  , 
V X y 


making  use  of  the  assumption  that  the  target  is  a body  of  rotation  about 
the  z-axis  (see  Appendix  8.2). 


a *(a  xa)=un  - nu 
u n X y z y z 


a • (a  X a ) = V n 
V n X y z 


a X a = (u  eos  a - u sin  a , - u eos  a , u sin  a) 
u r y z X X 
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a X 

V 


r 


o 


(v  cos  'X 

y 


-V  cos  a 

X 


V sin  a) 

X 


Consider  now  two  points  symmetric  about  the  yz  plane:  the  quantities  u , 

v^,  and  n^  have  opposite  signs.  Hence  neglecting  terms  that  will  cancel 
v;hen  integrating: 


line  u n - n u 

J ^ -,,,inc  /s  y z y z 

— ^ — X a = 2i!  a — ^ 

Y r X 

o 


Only  for  the  case  Cx  = 90  (broadside 


►me  ^ 

J X a ds 
r 
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2 

u cos  a-u  sin  a vn  cos  i 

^ 5 + _iL  z 

1 - 1 + C 


incidence) , this  reduces  to  simply 
2H^^^  a ds  proj  , 

X 

(1  - C) 


where 


ds  . = (a  • a ) dS  . 
prog  n 


I.n  the  above,  the  expressions  in  Appendix  8.2  were  used.  Using  a similar 

orocedure  we  find  that  for  the  other  oolarization  a = a : 

H p 


tine 


2H  a 


X a 


dS  = 


Ji 

(1  + ;) 


dS 


prog 


if  a:id  only  if  a = 90°. 

8.5  PROJECTED  AREA  OF  NEWLV  ILLUMINATED  REGION 

The  projected  area  of  tlio  nev/  region  needs  to  lie  related  to  its 
orientation  ia  order  to  determine  tlie  coordinates  and  derivatives  of  the 
new  region. 


The  region  is  illustrated  in  cross-section  and  perspective  in 
Fig.  41.  It  is  a slice  of  a skewband  about  the  body  of  rotation.  V.'e 
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r 


consider  first  only  one  side  of  the  new  region,  although  in  general  there 
will  be  two  disconnected  pieces. 

The  projected  height  at  a point  (y,z)  is  given  by 


h = 2/^-  ’ 

2 


where  is  the  radius  of  revolution  through  the  point  the  far  edge 

of  the  new  region.  Define  similarly  p^,  the  near  edge  of  the  region. 

Let  a be  the  surface  normal  at  x = 0 
n 

o 

a = (0,1,  -p')//!  + p'2' 
n 

o 


h(a  *a  ) ds  , 
n r 
o o 


where 

a = (0,  sin  cx,  cos  a)  is  the  direction  to  the  far  field, 
r 

o 


rhis  expression  needs  to  be  simplified. 


The  average  value  of  h is 


h = T (h,  -h_)  = /'(p,  + p')  (P,  -P,)'  ~ /2pAp'' 
av212  1212 

where  p is  the  average  of  P^'P2  Ap  = p^  - p^.  If  p is  taken  as  the 

value  in  the  region  center,  h is  then  an  underestimate  of  the  effective  h. 

av 

An  upper  limit  on  the  effective  h is 


h^  = 2/p^  - p^  = 2/pAp 
t 2 
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where  again  p is  the  value  in  the  center. 


Here  we  take  an  intermediate  value 


^eff  = 7 ' 


which  is  also  the  value  obtained  for  a sphere,  which  can  be  calculated  ex- 
actly. Also  note  that 


5s  = dz  /l  + p'^ 


Hence  obtain 


As  • = ■?■  /(^p'  (Az  sin  a + Ap  cos  a)  , 
pro3  2 

where  Ap  = p^  - p^  and  Az  = z^  - . 

In  the  case  of  general  oblique  incidence,  the  As  on  the  two  sides 
are  different,  but  in  the  case  of  a = 90°  we  multiply  the  above  by  2,  to 
account  for  both  sides.  For  the  case  of  cx=90°,  including  both  pieces  of 


the  new  region,  we  write 


i.e.  we  can  obtain  Az  and  hence  p^  = Ap/Az  since  known 

from  the  far  field  response. 


In  the  case  of  a = 0 , it  is  simpler  to  use  the  fact  of  rotational 


symmetry  and  write 


^Voj  = ^^^2  - 


P = -’^Pl  WPi  + 
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For  general  oblique  incidence,  the  main  problem  is  in  determining  the  rela- 
tive size  of  the  regions  on  each  side. 

8.6  INITIAL  RAMP  RESPONSE  OF  SPECULAR  POINT 

As  a first  step  in  the  inverse  scattering  process,  the  leading 
edge  of  the  ramp  response  is  analyzed  in  order  to  obtain  the  geometric  char- 
acteristics of  the  specular  point.  The  response  of  a curvilinear  patch  was 
derived  in  Section  5.3.  Specializing  those  results  to  the  specular  point 
and  writing  S for  the  projected  area  and  t for  the  time  since  incidence 
gives : 


r H (t) 
o r 


^ I t + — i — ± ~ — I — 

27T  dt  ^o  ^ 2tt  , 2 'o  2 " 47T  dt  'o  2 


/K 


dt 


(K  - K ) 
U V 

2 


(77) 


The  first  two  terms  are  the  projected  area  to  second  order  the  last  terms 
are  the  polarization  dependent  physical  optics  corrections,  where  the  sign 
is  + for  X polarization  and  - for  y polarization. 

The  derivation  of  the  above  assumed  that  both  curvatures  wore  non- 
zero. Additional  terms  arise  if  only  one  of  the  curvatures  is  non-zero  as 
in  the  case  of  a cylinder  viewed  at  broadside. 

The  projected  area  of  a cylinder  viewed  at  broadside  is 


S = 2L  / t (a  - t/4)'  , 

where  a is  the  radius.  At  t close  to  zero  we  can  approximate 

dS  = 2L  /a  t “ dt. 

Wo  use  this  expression  for  df>  in  Eq.  (48)  for  tlie  integral  over  tlie  self- 
patch  to  obtain  the  next,  [xslarization  dependent  term: 
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r H (t)  = — — t' 
or  '1 


~r  + r 

oTT»  a or 


C^H) 


1 CJJ’ 

In  the  above,  we  have  used  the  tacts  that  K =0,  r = — , ami  wheie  H js 

u V a i: 

tlu'  initial  rain[)  lesiJOtisi’  of  the  caps  as  qiven  iiy  Kq.  ('ll).  The  + siqn  in 
the  polarization  tiepi'iident  term  refers  to  tlie  case  of  x-poiar  i zat  ion , tin-  - 
siqn  to  p-polar izat ion.  ‘ 

In  tiie  inverse  piocedure  we  taki‘  tlie  averaqe^  of  t!ie  re.iponses  at 
the  two  polarizations  and  tost  for  tlio  presence  of  a dependence.  T)ie  co- 
efficii'nt  of  this  dei'cndenco,  by  Kq.  (7B)  , is  I'Xiiected  to  bo 


b >'  a 
71 


(79) 


Prom  the  remaininy  terms  due  to  tlie  caps,  the  two  principal  curvatures  at 
the  cap-side  of  the  boundary,  and  hence,  tlie  radius  a is  determined  by  the 
methods  of  Section  5.3.  P,q.  (79)  then  yives  the  lencjt.h  of  the  cylinder 
body . 


For  reference,  the  ideal  specular  response  of  a prolate  spheroid 
is  evaluated.  The  area  function  of  a prolate  spheroid  viev\ied  at  angle  i 
with  respect  to  the  z-axis  (axis  of  rotational  symmetry)  is 

3(t)  = Tia^blti'p  - tV4l/b^  ^ , 


v/here 


t = 0 at  time  of  incidence  at  the  specular  point, 
a,b  are  tiie  ellipse  axes 
b is  the  axis  of  rotation  (z) 

'2. 

H = b*'  cos  u + a^  sin'"  i . 

Taking  the  lierivatives  of  the  above,  evaluating  at  1 = 0 aiil  subs t i t ut  i na 
in  F.q.  (77)  yields 
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r H (t) 
o r 


2B 


a^b 


■K  - K > : 

U V ' t 


4B 


3/ 2 


2B 


; 2 


r,  '^K  - K \ 2 

1 , 1 i ^ + 

1 2 . 2 I 2B 


K -K  2^  ^3 

u V , a b \ t 


/'  4B 


3/2 


The  top  signs  here  refer  to  the  x-polarization. 


At  broadside  (a =90°),  which  is  the  test  case  of  this  report, 
above  reduces  to 


r H (t) 
o r 


^ + b 

4a  2 2 


_b 
4a  ' 


6 


For  a sphere,  both  polarizations  reduce  to 


r H (t) 
o r 


1 t 
4 


2 

2 ■ 


(80) 


the 


(81) 


(82) 
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MISSION 

of 

Rome  Air  Devebjrment  Center 


RADC  plans  and  conducts  research,  exploratory  and  advanced 
development  programs  in  command,  control,  and  communications 
(C^)  activities , and  in  the  areas  of  information  sciences 
and  intelligence.  The  principal  technical  mission  areas 
are  communications,  electromagnetic  guidance  and  control, 
surveillance  of  ground  and  aerospace  objects,  intelligence 
data  collection  and  handling,  information  system  technology, 
ionospheric  propagation,  solid  state  sciences,  microwave 
physics  and  electronic  reliability,  maintainability  and 
compatibility . 


